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DRSAR-PEL 


13  November  1979 


MEMORANDUM  FOR  RECORD 

SUBJECT:  Methods  for  Processing  Signals  from  Rotary  Potentiometers  with 
Application  to  the  Reduction  of  M110A2  Loader-Rammer  Test  Data 


1 .  Reference: 

a.  MFR,  DRSAR-PEL,  29  Aug  79,  subject:  Observations  Concerning  the 
M110A1  Loader-Rammer  Performance  Test  at  YPG,  Aug  79. 

b.  Article  appearing  in  the  Journal  of  Geophysical  Research,  Vol  77, 
No.  36  by  W.  D.  Hibler  III,  20  Dec  72,  title:  Removal  of  Aircraft  Alti¬ 
tude  Variation  from  Laser  Profiles  of  the  Arctic  Ice  Pack. 

2.  Background 

The  work  reported  in  this  memorandum  completes  one  task  in  a  project 
associated  with  a  problem  experienced  by  the  M110A2  SP  howitzer  system. 
Cannon  damage  reported  after  fielding  the  M110A2  in  Europe  in  1978  was 
found  to  be  associated  with  firing  the  M106  projectile  from  an  unseated 
position.  Some  projectiles  had  apparently  been  improperly  rammed  and  had 
fallen  back  upon  the  propelling  charge  when  the  gun  tube  was  elevated. 
Because  the  loader-rammer  (L/R)  was  implicated  in  this  problem,  a  series 
of  tests  were  devised  to  explore  the  limits  of  rammer  function.  Certain 
proposals  to  improve  the  standard  loader-rammer  were  submitted  and  as  a 
consequence,  the  scope  of  the  loader-rammer  tests  was  expanded  in  order 
to  compare  the  performance  of  the  standard  configuration  with  the  pro¬ 
posed  modification.  Tests  have  been  conducted  at  Yuma  Proving  Ground 
(YPG),  as  reported  in  Reference  a.  An  additional  L/R  test  program  was 
conducted  at  Jefferson  Proving  Ground  (JPG)  during  August  1979  using  a 
different  M110A2  weapon.  This  test  was  considerably  more  modest  in  scope 
than  the  Yuma  tests,  having  been  intended  as  a  corroborative  test.  The 
JPG  tests  measured  only  displacement  of  the  rammer  head  using  a  rotary 
potentiometer  (pot)  attached  to  the  axle  of  the  chain  sprocket.  From 
these  data  it  was  desired  to  derive  velocity  and  acceleration  of  the 
rammer  as  functions  of  time.  The  velocity  functions  for  different  experi¬ 
mental  systems  (or  treatments)  would  then  be  compared  to  identify  differ¬ 
ential  effects. 
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3.  In  practice,  however,  rotary  potentiometer  signals  contain  spurious 
noise  and  end-of-cycle  non-linearity.  Additionally,  the  displacement 
cycles  must  be  smoothly  pieced  together,  i.e.,  "stacked",  to  yield  a 
continuous  displacement  function  which  can  be  differentiated.  These 
effectsU)  (and  some  others)  pose  analytic  difficulties  which  must  be 
addressed  to  obtain  meaningful  velocity  functions. 

4.  Purpose 

There  are  two  principal  purposes  in  writing  this  MFR:  to  set  down 
the  lessons  learned  from  trying  various  methods  of  analysis  and  to  indi¬ 
cate  what  differences  appear  in  the  velocity  functions  estimated  for 
various  experimental  systems  run  at  JPG.  Relative  to  the  first  purpose, 
the  authors  would  be  gratified  if  some  of  the  techniques  used  in  the 
present  analysis  were  applied  or  refined  in  future  analyses  by  others 
confronted  with  similar  problems. 

5.  Organization 

The  remainder  of  this  memorandum  is  organized  as  follows:  Data  reduc 
tion  problems  and  methods  of  treatment  are  contained  under  Methodology. 
Results  of  applying  these  methods  to  the  present  problem  are  presented 
under  Resul ts ,  where  some  of  the  lessons  are  illustrated.  Finally,  the 
methodological  conclusions  and  specific  inferences  for  the  JPG  loader- 
rammer  data  are  presented  under  Conclusions.  For  the  reader  who  is  not 
interested  in  supporting  evidence  the  final  section  may  alone  serve  as 
an  adequate  summary. 

6.  Methodology 

The  displacement  of  the  rammer  head  as  measured  by  the  rotary  poten¬ 
tiometer  is  divided  into  cycles  each  of  which  represents  one  rotation  of 
the  pot.  The  potentiometer  signal  rises  to  a  maximum  voltage  and  then 
abruptly  falls  to  its  minimum  value  at  the  end  of  a  cycle.  However, 
experience  indicates  that  the  potentiometer  has  a  finite  resolution. 

This  causes  non-linearity  and  uncertainty  in  voltage  at  the  beginning 
and  end  of  each  cycle  and  produces  a  finite  fall  time  in  passing  from 


(1)  In  view  of  all  of  the  problems  encountered  in  obtaining  a  valid  velo 
city  function  from  rotary  pot  signals,  one  may  ask  why  not  use  an  alterna 
tive  type  of  instrument  such  as  a  tachometer?  In  the  present  case,  the 
authorsdid  not  have  a  choice  of  instruments  and  had  to  use  the  available 
data. 
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max  to  min  values.  Generally  also,  the  voltage  swing  corresponding  to 

one  rotation  of  the  pot  is  not  known  precisely  but  must  be  estimated  from 

the  noisy  output  signal.  To  obtain  a  displacement  signal  one  must  scan 

the  (analog)  record  for  max  and  min  values  --  y  and  y  .  --or  some 

max  min 

average  max  and  min  over  cycles  —  ym3v  and  y  .  —  and  obtain  a  scale 

Wax  mi n 

factor  of  displacement  per  volt  by  dividing  the  cycle  length  (24  inches 

in  the  present  case)  by  the  voltage  swing  (ym3v  -  y.  ).  This  scheme 

max  min 

is  implemented  in  the  computer  program  provided  in  Annex  2. 

7.  Processing  of  the  signals  will  generally  be  performed  on  a  digital 
computer,  so  that  an  analog-to-digital  (A  to  D)  conversion  of  the  pot 
signals  is  required.  Our  experience  indicates  that  this  operation  is 
also  a  source  of  noise.  In  fact,  A  to  D  conversion  seems  to  yield  quan¬ 
tization  of  the  signal  in  which  certain  "favorite"  numbers  repeatedly 
appear.  Ways  to  minimize  this  effect  include  recording  at  a  maximal 
signal  level  and  time  sampling  at  a  high  data  rate.  This  rate  may  be 
limited  by  the  available  high-speed  data  storage  but  should  be  at  least 
20  times  the  highest  significant*  frequency  seen  in  the  signal.  For 
example,  our  analog  pot  signals  were  time-sampled  at  500  hertz  (2  ms 
intervals)  since  an  upper  frequency  limit  of  25  hertz  was  expected.  Pre¬ 
sent  experience  indicates  that  filtering  the  analog  signal  with  a  100 
hertz  pass  band  before  digitizing  creates  more  problems  than  it  solves. 

This  type  of  filtering  exaggerates  the  non-linearity  in  the  end-of-cycle 
pot  signal  and  leads  to  stacking  difficulty. 

8.  To  obtain  a  stacked  displacement  signal  which  is  sufficiently  smooth 
to  tolerate  a  differencing  procedure  requires  some  digital  filtering. 

Two  types  of  non-recursive  filters  were  tried  —  a  symmetric  moving 
average  and  a  filter  having  exceptional  discrimination  and  low  ripple 
outside  the  passband.  The  latter  was  designed  by  a  method  devised  by 

W.  D.  Hibler  (Reference  b).  The  modulus  and  squared  modulus  of  the 
transfer  function  of  these  filters  are  shown  in  Figures  1  through  4  of 
Annex  1.  In  spite  of  the  theoretical  advantage  of  Hibler ‘s  filter,  it 
did  not  prove  to  be  as  satisfactory  in  this  application  as  a  simple  sym¬ 
metric  moving  average.  Consequently,  the  results  displayed  in  the  figures 
of  Annex  1  were  generated  using  a  symmetric  moving  average  filter  having 
a  maximum  absolute  lag  of  25.  With  a  sampling  period  of  2  ms,  this  moving- 
average  filter  has  a  resolution  bandwidth  of  10  hz  (equivalently,  an  aver¬ 
aging  period  of  0.1  sec).  An  averaging  period  of  about  0.1  sec  is  required 


*  At  the  -20  db  level . 
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to  achieve  the  proper  degree  of  smoothing  of  the  displacement  judged  nec¬ 
essary  for  this  application. 

9.  When  successive  potentiometer  cycles  are  stacked  by  adding  the  pre¬ 
vious  cumulative  full-cycle  value  to  the  pot  reading, smal 1  discontinuities 
appear  at  the  ends  of  cycles.  To  eliminate  these  and  other  outliers  from 
the  displacement  record  it  was  desirable  to  use  an  outlier  detecting  and 
purging  algorithm.  Simply  relying  on  filtering  to  reduce  the  effect  of 
these  outliers  was  not  practical.  The  outlier-detecting  procedure  works 
as  follows:  First,  the  unmodified,  stacked  displacement  signal  is  fitted 
with  a  (6th  degree)  polynomial  function  of  time  using  multiple  linear 
regression.  The  standard  error  of  the  estimate  from  this  regression  is 
then  used  as  a  measure  of  discrepancy  of  a  data  point  from  the  trend. 

Data  points  which  are  more  than  about  three  standard  errors  from  the  fitted 
function  were  replaced  by  the  value  of  the  fitted  function.  Finally,  the 
moving  average  was  calculated  to  attenuate  the  pot  noise.  Actually, 
expurgation  of  outliers  and  filtering  can  occur  in  the  same  DO-loop  of  a 
computer  program  as  was  done  in  the  program  in  Annex  2. 

10.  Having  obtained  a  suitably  smooth  displacement  function,  velocity  is 
estimated  using  a  first-order  central -difference  approximation: 

dy^/dt  =  (yi+]  -  yi_] )/2h, 

with  time  step  h. 

Although  the  values  of  y.  above  may  be  moving-average  displacements,  the 

differencing  operation  produces  high-frequency  noise  which  may  be  distract¬ 
ing  to  the  person  who  examines  this  unsmoothed  velocity  estimate.  Conse¬ 
quently,  we  have  found  it  convenient  to  smooth  the  above  estimate  by  pass¬ 
ing  the  unsmoothed  velocity  through  a  digital  filter  whose  passband  some¬ 
what  exceeds  that  of  the  moving-average  filter  which  was  applied  to  the 
displacement  signal.  In  the  present  application,  a  non-recursive  filter 
with  a  max  lag  of  6  was  used  for  this  purpose. 

11.  Another,  completely  different  approach  to  estimating  velocity  and 
acceleration  was  used  and  found  to  produce  valid  results  if  properly 
restricted.  This  approach  is  referred  to  as  analytic  since  the  velocity 
and  acceleration  can  be  represented  analytically  as  polynomial  functions 
of  time.  In  our  algorithm  the  movi ng-average  displacement  is  fitted  with 

a  sixth  degree  polynomial.  (A  fifth  degree  polynomial  works  almost  as  well.) 
This  function  is,  of  course,  analytically  differentiable  yielding  a  fourth 
degree  polynomial  estimate  of  acceleration.  Both  velocity  and  acceleration 
estimates  which  are  derived  in  this  manner  are  displayed  among  the  results 
in  Annex  1.  Clearly,  analytic  estimates  may  be  rather  gross  if  the  domain 
of  the  fitted  function  is  too  large  or  if  time  derivatives  beyond  the  sixth 
are  numerically  significant  anywhere  over  the  domain. 
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12.  Results 


The  results  of  this  study  are  presented  in  the  form  of  plots  of  out¬ 
put  from  the  computer  program.  The  same  types  of  graphs  are  given  for 
each  of  24  rammer  tests.  The  sequence  of  the  plots  is  as  follows: 

(a)  Raw  potentiometer  signal  (a  function  of  time); 

(b)  Displacement  signal  produced  by  stacking  the  pot  signal,  purging 
and  replacing  outliers  and  scaling; 

(c)  Analytic  estimate  of  velocity  using  the  6th  degree  polynomial  fit 
to  the  displacement  data  in  (b); 

(d)  Analytic  estimate  of  acceleration; 

(e)  Unsmoothed  estimate  of  velocity  obtained  by  differencing  the 
moving  average  displacement; 

(f)  Smoothed  velocity  estimate;  from  (e); 

(g)  Crossplot  of  the  smoothed  velocity  versus  the  analytic  estimate  of 
displacement. 

The  last  estimate  was  the  most  useful  in  comparing  the  results  from  different 
rams  since  it  is  independent  of  time  zero  which  generally  varied  between  rams 
The  sequence  of  graphs  proceeds  from  the  bottom  to  the  tO£  of  each  page 
through  all  ram  numbers  without  a  page  break. 

13.  A  description  of  the  experimental  systems  tested  is  provided  in  Table 
1.  Table  1  also  summarizes  the  average  (AVG)  and  standard  deviation  (SD) 
of  the  ram  speed  obtained  in  the  last  full  cycle  of  rotary  pot  motion.  The 
first  10  rams  were  conducted  using  the  M2A2  cannon,  i.e.,  an  M110  howitzer 
configuration.  The  remaining  14  rams  were  conducted  in  the  M110A2  howitzer 
using  an  M201  cannon.  Because  of  the  difference  in  cannons  the  total  travel 
of  the  rammer  head  is  85  inches  in  the  MHO  and  about  93  inches  in  the  M110A2 
Due  to  this  difference  the  same  rammer  produces  a  somewhat  differently  shaped 
velocity  function  in  the  two  systems.  This  difference  is  quantified  in  the 
Cone! usions .  The  standard  loader-rammer  (L/R  STD)  was  used  for  the  first 
five  rams  and  the  modified  configuration  (L/R  MOD)  for  the  next  five.  As 
noted  in  Table  1,  both-  L/R  configurations  were  also  tested  with  the  M201 
cannon.  For  each  of  these  configurations  testing  was  done  at  maximum 
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acceptable  L/R  system  pressure  (5  rams  each)  and  at  minimum  acceptable 
L/R  system  pressure  (2  rams  each).  Although  velocity  differences  between 
L/R  configurations  in  the  M110A2  are  slight,  the  effect  of  L/R  system 
pressure  on  velocity  is  noticeable.  A  reduction  in  the  ram  speed  near 
end-of-travel  of  from  0.4  to  1.0  f/s  accompanies  the  noted  pressure  reduction. 

14.  The  data  set  obtained  from  digitizing  the  analog  rotary  pot  signal  for 
a  typical  ram  consists  of  approximately  600  points.  Sets  of  this  size  per¬ 
mit  data  processing  of  the  sort  described  above  while  retaining  an  adequately 
large  value  of  degrees  of  freedom.  Of  course,  the  use  of  a  non-recursive 
filter  shortens  the  unfiltered  data  set  at  each  end  by  the  value  of  the 
maximum  lag  used  in  the  filter.  Thus,  in  this  application  25  data  points 
were  eliminated  from  each  end  of  the  displacement  data  to  obtain  the  moving 
average  displacement. 

15.  The  first  twenty  figures  in  Annex  1  are  intended  to  illustrate  the 
effect  of  various  errors  associated  with  the  data  and  with  the  processing 
methods.  These  illustrations  support  previous  claims  of  lessons  learned 
and  give  the  reader  some  "feel"  for  the  problems  encountered  in  using  data 
from  rotary  potentiometers.  Starting  with  Rammer  Test  No  1  shown  in  Fig 

5,  one  can  observe  the  kind  of  stacking  error  typically  encountered  in  pro¬ 
cessing  the  pot  signal.  Without  identifying  and  purging  the  resulting  out¬ 
liers  from  the  displacement  signal  (Fig.  5)  a  characteristic  set  of  spikes 
is  produced  in  the  unsmoothed  velocity  estimate  (Fig.  7).  Even  after  smooth¬ 
ing  the  velocity,  anomalous  bumps  (reduced  spikes)  are  seen  in  the  velocity¬ 
time  plot  (Fig.  8)  and  velocity-displacement  plot  (Fig.  9).  One  can  also 
see  the  effect  on  velocity  of  potentiometer  non-linearity  --  the  0.1  sec 
shallow  waves  appearing  in  the  velocity  after  level  off. 

16.  Figure  10  illustrates  another  approach  to  filtering  the  displacement 
signal  prior  to  differencing  to  obtain  velocity.  The  Hibler  filter,  describ¬ 
ed  earlier,  was  applied  to  the  unexpurgated  displacement  (Fig.  5)  and  differ¬ 
enced  to  yield  the  unfiltered  velocity  shown  in  Figure  10.  Because  of  the 
sharp  frequency-domain  cutoff  and  low  ripple  in  the  stop  band,  one  would 
expect  the  Hibler  filter  to  perform  better  than  a  simple  moving  average. 
However,  the  smoothed  velocity  estimate  (Fig.  11)  is  not  subjectively  better 
than  its  moving-average  counterpart  (Fig.  9).  Further  reduction  of  the 
cutoff  frequency  of  the  filter  would  improve  smoothness  and  reduce  the 
effect  of  spurious  spikes  in  the  displacement  but  would  lack  resolution  and 
introduce  velocity  bias.  A  better  strategy  for  data  reduction  is  to  prepro¬ 
cess  the  unfiltered  displacement  to  detect  and  purge  outliers.  This  proce¬ 
dure  does  not  sacrifice  resolution  bandwidth  as  does  narrow-passband  filter¬ 
ing.  See  Figure  19  for  a  better  estimate  using  data  from  ram  test  number  1. 
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Figure  12  also  illustrates  the  effect  of  smoothing  with  Hibler's  filter  -- 
in  this  case  on  the  data  from  ram  number  2.  Here  also,  one  may  conclude 
that  the  velocity  estimated  from  differencing  the  filtered  displacement 
is  not  better  than  simply  using  a  moving  average  on  the  expurgated  data. 

See  Figure  29  for  comparison. 

17.  A  surprising  thing  happened  to  the  data  for  several  rams:  An  A  to 
D  transcription  and/or  tape  copying  error  caused  intervals  of  data  to  be 
omitted  in  the  final  (digital)  record.  Figure  13  illustrates  this  incident 
in  the  pot  signal  for  ram  test  number  5  and  Figure  14  indicates  the  result¬ 
ing  error  in  the  expurgated  displacement.  Evidently  the  outlier-detecting 
algorithm  cannot  cope  with  this  kind  of  error.  The  consequence  of  this 
error  is  serious  as  can  be  seen  in  Figures  15  and  16.  In  this  case  some 
manual  editing  of  the  data  was  necessary  before  reprocessing  (Figs.  43  and 
44). 


18.  One  additional  type  of  error  is  noteworthy  --  inclusion  of  too  many 
pre-ram  and  post-ram  data  points  in  the  displacement  record.  If  a  long 
post-ram  plateau  occurs,  as  in  Figure  23,  the  analytic  velocity  estimate 
(Fig.  24)  will  be  quite  poor.  The  sixth-degree-polynomial  fit  becomes 
inadequate  if  the  domain  is  not  properly  bounded.  By  contrast,  a  displace¬ 
ment  record  having  a  truncated  plateau  (Fig.  25)  yields  an  excellent  velo¬ 
city  estimate  (Fig.  26)  and  a  reasonable  acceleration  function  (Fig.  27). 

No  attempt  was  made  to  rectify  the  analytic  estimate  of  acceleration  out¬ 
side  of  its  applicable  domain. 

19.  One  is  able  to  observe  specific  differences  in  the  ram  velocity  func¬ 
tions  between  different  experimental  systems  (treatments)  even  when  the 
data  contain  the  imperfections  described  here.  One  reason  for  this  discrim¬ 
ination  is  the  consistent  manner  with  which  stacking  error  and  filtering 
bias  enter  the  velocity  estimated  from  run  to  run.  These  velocity  differ¬ 
ences  in  the  JPG  data  are  discussed  under  Conclusions. 

20.  Conclusions 


Conclusions  are  summarized  here  under  two  categories:  general  metho¬ 
dological  observations  and  specific  inferences  about  the  Ml  10/Ml 10A2  L/R 
tests  at  JPG.  Relative  to  methods  for  analyzing  rotary  potentiometer 
signals: 

(1)  It  is  best  to  stack  digitized  rotary  potentiometer  signals  which 
have  an  abrupt  drop  from  max  to  min  values.  Therefore,  a  broad-band  analog 
signal  is  a  desirable  starting  point  since  this  affords  the  most  abrupt  drop. 
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(2)  Some  stacking  error  appears  to  be  unavoidable.  Therefore,  it  is 
essential  to  process  the  stacked  displacement  signal  with  an  outl ier-detect 
ing  and  purging  algorithm.  As  a  minimum,  such  an  algorithm  should  obtain 
departures  from  a  global  average,  as  was  done  here  by  polynomial  regression 

(3)  Occasional  analog-to-digital  encoding  errors  and  digital  data 
omissions  occur.  Therefore,  it  is  important  to  display  the  digitized 
potentiometer  signal  graphically  to  assist  in  detecting  this  type  of 
error.  In  this  case  manual  editing  of  the  data  is  necessary. 

(4)  Elimination  of  stacking  and  encoding  errors  solely  by  high- 
discrimination  digital  filters  (such  as  Hibler's)  seems  to  be  impractical, 
since  a  sufficiently  narrow  filter  frequency  passband  (or  long  period)  to 
reduce  the  effect  of  the  errors  in  the  output  signal  would  incur  an  objec¬ 
tionable  degree  of  bias. 

(5)  Although  somewhat  cosmetic,  applying  a  digital  filter  to  the 
velocity  signal  to  smooth  it  is  recommended.  To  avoid  additional  bias 
the  passband  of  this  filter  should  be  larger  than  that  applied  to  the 
displacement  signal  used  to  develop  the  velocity  estimate. 

(6)  If  a  gross-average  velocity  estimate  is  desired,  one  can  obtain 

a  good  analytic  approximation  from  a  high-degree  polynomial  in  time  fitted 
to  the  displacement  data  over  a  restricted  domain  of  this  function.  Do 
not  expect  this  analytic  velocity  estimate  to  hold  outside  of  this  limited 
domain. 

(7)  Because  of  the  generally  arbitrary  nature  of  time  zero  and  be¬ 
cause  of  temporal  shifts  due  to  filtering,  it  is  recommended  that  compari¬ 
sons  of  two  distinct  runs  be  made  in  phase  space,  i.e.,  via  a  crossplot 

of  velocity  vs  displacement. 

(8)  In  making  phase-space  plots  it  is  desirable  to  use  a  highly 
smoothed  estimate  of  displacement,  e.g.,  by  using  regression,  so  as  not 
to  incur  anomalous  multi-point  function  incidents. 

21.  With  respect  to  the  results  of  the  Mil  0/Ml 10A2  L/R  tests  at  JPG: 

(1)  The  form  of  the  velocity  versus  displacement  plots  are  nearly 
identical  for  repeated  rams  with  a  specific  experimental  system.—  cannon 
type,  L/R  type,  and  L/R  system  pressure. 

(2)  Differences  in  these  plots  appear  between  experimental  systems  in 
two  respects:  the  max  velocity  level  achieved  and  the  displacement  at 
which  a  given  level  is  achieved. 
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(3)  During  the  first  ten  ramming  tests  using  the  M2A2  cannon  at  max 
L/R  pressure,  rams  1  through  5  apply  to  the  L/R  STD  configuration  and  rams 

6  through  10  to  the  L/R  MOD.  The  former  generally  achieve  a  slightly  higher 
peak  velocity  than  the  latter.  Another  distinction  in  these  tests  is  in 
the  shape  of  the  velocity  function.  The  velocity  of  L/R  STD  exhibits  a 
somewhat  more  convex  (hump-backed)  shape  than  that  of  the  L/R  MOD. 

(4)  The  shape  of  the  velocity  function  which  characterizes  rams  1 
through  5  also  differs  from  that  of  rams  11  through  15.  The  latter  set 
has  a  decidedly  sway-back  shape  with  a  shallow  dip  in  the  velocity  at 
midram.  Additionally,  a  somewhat  shorter  rise  to  peak  velocity  is  noted 
in  the  latter  set.  The  difference  in  experimental  system  in  these  sets  is 
simply  a  difference  in  cannon  type;  the  M2A2  cannon  was  used  in  rams  1 
through  5  whereas  the  M201  cannon  was  used  in  rams  11  through  15.  Apparent¬ 
ly,  the  standard  rammer  tends  to  increase  its  speed  slightly  at  the  longer 
travel  experienced  in  the  M201  cannon  --  93  inches  versus  85  inches  in  the 
M2A2  cannon. 

.  (5)  A  reduction  in  the  L/R  system  pressure  from  max  to  min  acceptable 
under  otherwise  identical  conditions  causes  a  reduction  in  peak  ramming 
speed.  This  reduction  is  anticipated  because  the  level -off  speed  is  con¬ 
trolled  by  a  restricted  flow  of  oil  under  a  pressure  differential  which  is 
nearly  constant  throughout  a  ram  cycle.  In  the  L/R  STD  configuration  this 
speed  reduction  amounted  to  only  about  0.4  f/s  based  upon  the  average  over 
the  last  complete  pot  cycle  and  about  1  f/s  based  upon  the  difference  in 
maxima  (rams  11  through  15  vs  16  and  17)  whereas  in  the  L/R  MOD  configura¬ 
tion  this  reduction  was  about  1  f/s  (rams  18  through  22  vs  23  and  24). 

For  both  configurations  the  shape  of  the  velocity  curve  in  phase  space  was 
not  altered  by  the  reduction  in  L/R  system  pressure. 


2  Incl 

1.  Annex  1 

2.  Annex  2 


GEORGE  SCHLENKER 
Operations  Research  Analyst 


LANNY  WELls 

Operations  Research  Analyst 


15 


' 


16 


ANNEX  1 


RESULTS  FROM  M110/M110A2 
LOADER-RAMMER  TESTS  AT 
JEFFERSON  PROVING  GROUND, 
AUGUST  1979 
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TABLE  1.  AVERAGE  FINAL ^ '  RAMMING  SPEED 
FOR  THE  M110A2  LOADER-RAMMER 
OBTAINED  IN  THE  JPG  TESTS,  11  AUG  79 


ROUND 

NUMBER 

SPEED 

t£ZSi 

COMMENTS 

1 

9.98 

M2A2  Cannon 

2 

9.90 

L/R  STD(2) 

3 

9.80 

System  Pressure 

4 

9.76 

5 

9.85 

AVG/SD 

9.858/0.086 

6 

9.55 

M2A2  Cannon 

7 

9.44 

L/R  MOD ( 2 ) 

8 

9.44 

System  Pressure 

9 

9.41 

10 

9.45 

AVG/SD 

9.458/0.054 

11 

9.70 

M201  Cannon 

12 

9.59 

L/R  STD' 2) 

13 

9.58 

System  Pressure 

14 

9.66 

15 

9.44 

AVG/SD 

9.594/0.099 

16 

9.50 

System  Pressure 

17 

9.28 

AVG/SD 

9.39/0.156 

18 

9.68 

M201  Cannon 

19 

9.89 

L/R  MODI* 1 2 3) 

20 

9.80 

System  Pressure 

21 

9.79 

22 

9.67 

AVG/SD 

9.766/0.092 

23 

8.65 

System  Pressure 

24 

8.84 

AVG/SD 

8.745/0.134 

Max^3^ 


Max^ 


Max^3^ 


Max^3^ 
Mi  n  ^ 3  ^ 


Notes: 

(1)  The  "final"  ramming  speed  was  calculated  as  the  average  over  the 
last  complete  rotary  pot  cycle. 

(2)  Two  loader-rammer  configurations  were  tested  --  the  current  stan¬ 
dard  (L/R  STD)  and  a  proposed  modification  (L/R  MOD). 

(3)  The  pressure  in  the  oil  or  nitrogen  of  the  loader-rammer  system 
was  adjusted  prior  to  each  ram  using  the  gage  provided  with  the 
system. 
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FIGURE  1 .  MODULUS  OF  THE  TRANSFER  FUNCTION  OF  A 
SYMMETRIC  MOVING  AVERAGE  DIGITAL  FILTER 
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FIGURE  2.  SQUARED  MODULUS  OF  THE  TRANSFER 
FUNCTION  OF  A  SYMMETRIC  MOVING 
AVERAGE  DIGITAL  FILTER 
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FIGURE  3.  MODULUS  OF  THE  TRANSFER  FUNCTION 
OF  A  SYMMETRIC  NONRECURSIVE  DIGITAL 
FILTER  USING  HIBLER'S  METHOD 
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FIGURE  4.  SQUARED  MODULUS  OF  THE  TRANSFER 
SO  MOOULUS  FUNCTION  OF  A  SYMMETRIC  NONRECURSIVE 
O  DIGITAL  FILTER  USING  HIBLER'S  METHOO 
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RAMMER  TEST  01  AT  J.P.G.  FIGURE  6. 
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RAM  EP  TE5T  01  AT  j.P.G.  FIGURE  7. 


Velocity  estimate  obtained  from  the 
moving-average  displacement  record 
having  the  stacking  error.  Anomolous 
spikes  appear  at  the  ends  of  the  moving 
average  window. 
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RAMMER  TEST  01  AT  J.P.G.  FIGURE  9. 


RAMMER  TEST  01  AT  J.P.G.  FIGURE  8. 
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RAMMER  TEST  01  AT  J.P.Q.  FIGURE  11. 
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RAMMER  TEST  01  AT  J.P.G.  FIGURE  10. 
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RAMMER  TEST  02  AT  J.P.G.  FIGURE  12. 
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RAMMER  TEST  05  AT  J.P.G.  FIGURE  13. 
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RAMMER  TEST  05  AT  J.P.C. 


FIGURE  16. 
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RAMMER  TEST  05  AT  J.P.G.  FIGURE  15. 
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Illustration  of  the 
effect  of  missing  data 
on  the  velocity  estimated 
from  a  moving  average 
displacement  signal. 


^.00 


0.20 


0.60 


1  .00 


TIME  (SEC) 


O'.  40 


30 


0.8G 


o i spl  sig  mm  oispl  sic  mv) 


o 

o 


RAMMER  TEST  01  AT  J.P.O.  FIGURE  18. 
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RAMMER  TEST  01  AT  J.P.O.  FIGURE  17. 


TIME  ( SEC i 


31 


VEl.  (F/S)  VEl  (F/S) 


o 

o 

CM 


RAMMER  TEST  Oi  AT  J .P.Q.  FIGURE  20. 
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RAMMER  TEST  01  AT  J.P.G.  FIGURE  19. 
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RAMMER  TEST  02  AT  J.P.G.  FIGURE  22. 
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RAMMER  TEST  01  AT  J.P.G.  FIGURE  21 
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purged  potentiometer  data. 
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RAMMER  rE5T  02  AT  J.P.G.  FIGURE  24. 
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RAMMER  TEST  02  AT  J.P.G.  FIGURE  23. 
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RAMMER  TEST  02  AT  J.P.G.  FIGURE  26. 
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RAMMER  TEST  02  AT  J.P.G.  FIGURE  28. 


RAMMER  TEST  02  AT  J.P.G.  FIGURE  27. 
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S  RAMMER  TEST  02  AT  J.P.G.  FIGURE  30. 
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RAMMER  TEST  03  AT  J.P.G.  FIGURE  31. 
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FIGURE  34. 
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RAMMER  TEST  03  AT  J.P.Q.  FIGURE  33. 
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RAMMER  TEST  03  AT  J.P.O.  FIGURE  35. 
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RAMMER  TEST  04  AT  J.P.G.  FIGURE  38. 
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RfittflER  TEST  04  AT  J.P.G*  FIGURE  37. 
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ANNEX  2 

COMPUTER  PROGRAM  SOURCE 
LISTING  FOR  EST.VEL.ACC 

A  PROGRAM  FOR  PROCESSING 
DISPLACEMENT  SIGNALS  FROM 
ROTARY  POTENTIOMETERS  TO 
ESTIMATE  VELOCITY 
AND  ACCELERATION 


43 


44 


c 

C  PROGRAM  FOR  PROCESSING  DISPLACEMENT  SIGNALS  FROM  ROTARY 
C  POTENTIOMETERS  TO  ESTIMATE  VELOCITY*  AND  ACCELERATION 
C 

C  THIS  PROGRAM  CONTAINS  EXTERNAL  REFERENCES  TO  PRIME  400  COMPUTER 
C  SYSTEM  ROUTINES  AND  CRLCOMP  925/936  SOFTWARE-  ON  THE  ARRCOM  S4E 

C  SYSTEM  THE  FOLLONING  LIBRARIES  MUST  BE  LOAOEO: 

C  1.  LWVLIB 

C  2-  VAPPLB 

C  3-  VPLOTLB 

C  4.  VBPOOt 

C 

REAt*8  OATAC  1200/7)  ,X(7)  .FJ*.FJJM1 . VEL .ACC .S .RSQ .OISPL . VELT .OELT . 
♦RESIO.XT .YEST0.AC66) 

DIMENSION  Y ( 1200) ,YBAR( 1200) ,S2( 1200).INF0(  12)  ,N(  2)  .TM  1200) . 

♦VMAC 1200) rVMASC 1200)  .C0F2C 13) ,YEST( 1200)  JIPOTC 1200) 

INTEGER  TNC 16) .TITLE! 40) .COMOFLC 16) ,BUF( 60 ) 

LOGICAL  NSCALE 
C 

COMMON/CPLCOM/  AXL.AYL 
C 

C  THE  OATAC I. J)  ARRAY  IS  USEO  TO  STORE  THE  I 'TH  VALUE  OF  TIME  TO  THE 

C  J’TH  POWER  FOR  EXPONENTS  UP  TO  ANO  INCLUOING  NOEG.  THE  VECTOR 

C  NITH  JsNOEG+1  CONTAINS  THE  VALUES  OF  THE  SCRLEO  ANO  STRCKEO  OEPEN- 

C  OENT  VARIABLE.  DISPLACEMENT  (IN).  THE  UNSCALEO  ROTARY  POTENTIO- 

C  METER  ORTA  ARE  STOREO  IN  RPOT, 

C 

$  INSERT  SYSCOM>A*KEYS 
♦  INSERT  SYSCOM>KEYS >F 
♦INSERT  SYSCOM>ERRO.f 
C 

OATA  0ATA/8400*0. 0/ . INFO/1 ,2.  :100000.3*  ’  *  -6*0/ ,Vt1A/1200*0 .0/ , 

♦YE ST/ 1200* 0-0/ 

OATA  C0F2/1 .  1286702E-02 -3  -0793781E-02 .5 -7062223E-02 .8 -6293384E-02 . 

1  1 .1312581E-01 ,1 .3202178E-01 -1 -3883263E-01 .1 .3202178E-01 . 

2  1 .1312581E-01 .8.6293384E-02.5.7062223E-02.3.0793781E-02. 

3  1 .1286702E-02/ 

C 

C  CALL- OPEN^LC ’INPUT  FILE  TREENAME? . » ->  ’  .26 .A^REAO ,TN .32  » 1 ) 

C 

C  PROCESS  ORTA  DIRECTLY  FROM  TAPE 
C 
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C  CALL  TNOUAC  'TAPE  UNIT»?...>  ’.16) 

C  REAOC1.*)  NO 

NIL=4 

C  CALL  C*M13( l.O.NU.O)  /*  OPEN-  TAPE  TO  READ 

CALL  OPENtAC  AHREAO  ,  ’PEXLOH>A>TAPE-OUMP  M8  .NU) 

C 

00  104  NFILE=1 ,24 
C  CALL  C$M13( 6.0.NU .0 ) 

CALL  CHA007CNU) 

C  IFCNFILE.LT  .6 )  GO  TO  104 

C 

C  OPEN  SCRATCH  FILE  TO  SAVE  TERMINAL  OALOGUE 

C 

CALL  OPNSFLC  3  ,COMOFL  ) 

CALL  CLOS$A(  3 ) 

CALL  COMO$*(  :  000020 .COMOFL . 32 .0 .COOE ) 

C 

IF(COOE.GT.O)  GO  TO  103 

C  CALL  OPENtLC 'OUTPUT  FILE  TREENAME? .  .  .>  ’ .26 . A$NRIT ,TN , 32  ,2 ) 
ENCOOEC  32 . 12 . TN )  NFILE 
12  FORMATC  'RAMMER— TEST#  ’  ,B ’««»  * ) 

CALL  OPENtAC AtWRIT .TN.32.2) 

CALL.  TRNCtAC  2 )  /*  TRUNCATE  OUTPILT  FILE 

C 

C  READ  TITLE  OF  OATA  SET 

C 

C  REA0C5.2)  TITLE 

2  FORMATC  40A2 ) 

C  CALL  IHAM13CNU. TITLE. 40.0) 

CALL  I*A0D7(NU. TITLE. 40,0) 

NRITEC  6.1)  TITLE 
1  FORMAT (  1H  .40A2 ) 


C  CALL  I*AM13(NU. TITLE. 40.0) 
CALL  IHA0D7CNU. TITLE. 40.0) 
NRITEC  6.1)  TITLE 


c 

CALL 

CALL 

C$M13( 5 .0 . NU.O ) 

I$A0D7( NU .BUF . 40 .0 ) 

/* 

SKIP 

3R0  REC 

c 

CALL 

CALL 

lit  AM13CNU.&UF. 60.0) 

I $A007( NU.BUF.60.0) 

/* 

READ 

RECORO  INTO  BUF 

OECOOEC 120.11  ,BUF )  POINTS .OELTAT 
11  F0RMATC2G16.8) 

TIME0=0 .0 


/*  SKIP  TO  NEXT  FILE 

/*  SKIP  FIRST  4  FILES 
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NLZERO-D 

NDATA=PO I  MTS-1  -0 
NSET=NDATA 
YSCRLE=1D.9D9 
YMIN=-D .01 

IHIB=D  /* 

FC=5.0  /* 

YD=0.0 

NCH=NLEM$A( TITLE .4D ) 

C 

C  READ  COMTROL  P ARAMS 

C 

C  READC5.*)  TIMED. NLZERO .NSET .NDATA .DELTAT . YSCALE . YMIN . YO 
C  NSCALE=YSNO$A(  ’SELF  SCALING  DESIRED?... >  ’  .26.A$NDEF) 
NSCALE=. FALSE. 

DELTrDELTAT 
ISTART=NLZER0+1 
I STOP=NLZERO+NOATA 
C 

C  READ  THE  OUTLIER  LAG  PRRAM 

C 

C  CALL  TNOUAC ’OUTLIER  LAG  PARAMETER? .. •>  ’.27) 

C  READC1.*)  MLAG 

MLAG=25 
C 

C  READ  MIN  NUMBER  OF  STANDARO  DEVIATIONS  TO  DECLARE  AK  OUTLIER 
C 

C  CALL  TNOUAC ’NUMBER  OF  SO  S  FOR  OUTLIER?... >  ’.32) 

C  READ! 1 .» )  SDNO 

SDN0=2.64 
C 

C  READ  DEGREE  OF  POLYNOMIAL  FIT 

C 

C  CALL  TNOUAC  ’DEGREE  OF  FITTING  POLYNOMIAL? .. .>  ’.34) 

C  READ! 1 ►» )  NDEG 

NDEG=6 

NDEGP1=NDEG+1 
MAVG=2*MLAG+1 
FMAVG=FLOAT(  MAVG) 

LSTART=ISTART+MLAG 
LSTOP=ISTOP-MLAG 
I 1=ISTART+MRVG-1 


INDICATOR  VARIABLE  FOR  HIBLER’S  FILTER 
CUTOFF  FREQUENCY  (HZ)  IN  FILTER 


I _ • 
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C  REAOCS.*)  (  Y(  I ) .  1= I  START  « I STOP ) 

NREC=NDATA/8  +  C1IN0(  1  ,MOOC  NDflTfi  ,8 ) ) 

DO  13  I=1.NREC 

C  CALL  I$AM13( NU  »BUF .60 ►♦IB ) 

CALL  I$A007(MJ.BUF.60,$15) 

DO  TO  13 
15  CALL  PRERR 
CALL  EXIT 

C  IFC  NFILE  *NE  .29 )  DO  TO  105 
C  CALL  DETERRC  BUF .2 ) 

C  IFC BUEC 1  ).NE. ’IE’)  DO  TO  106 

C  CALL  I$AM13(  NU  .BUF  r60 .0 ) 

C  CALL  I $R007( NU  »BUF .60.0) 

13  DECODE ( 1 20 .* .BUF )  C Y( 8*I+ISTART-9+ J) . J=1 .8 ) 

C 

C  DETERMINE  AN  AVERAOE  INITIAL  VALUE 

C 

SUM=0.0 

ISP19=ISTART+19 
DO  94  I=ISTART .  ISP19 
94  SUM=SUM+Y(I) 

Y0=SUM/20 .0 
YLIM=93  .O/YSCALE+YO 
C 

YRANOE=24 .O/YSCALE  /*  24  INCHES  PER  POT  CYCLE 

YMAX=YMIU+YRANGE 
T0LER=0 .OlwYRANGE 
IFC  *NOT .NSCALE )  DO  TO  90 
C 

C  SCAN  DATA  FOR  MAX  AND  MIN 
C 

NCYCLE=1 
YMAX=YC ISTART) 

YMIN=YC ISTART) 

ID0=ISTART+1 
DO  80  1  =  100 .ISTOP 
IFCY(I).DT.YMAX)  YMAX=YCI) 

IFCY(I).LT.YMIN)  YMIN=YCI) 

80  CONTINUE 

YRANGE=YMAX-YMIN 
YSCALE=24 .O/YRANDE 
TEST=YRANGE/2.0 
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YIM1=Y( ISTRRT ) 

DO  82  1=100. ISTOP 

IF(R6S(YIM1-Y( I)) .OT.TEST)  NCYCLE=NCYCLE+1 
YIM1=Y( I ) 

82  CONTINUE 

DISPF=YSCALE*( Y(  ISTOP)-YD+Ft0AT(  NCYCLE-1  )*YRRNGE  ) 
IF(DISPF.0T.83.D  ^ND.  DISPF  .LT  .94 .0)  GO  TO  9D 
C 

C  FIRST  SELF-SCALINO  PROCEDURE  PRODUCES  FINAL  DISPLACEMENT  OOT-OF- 
C  BOUNDS.  TRY  A  SECOND. 

C 

NRITE( 1 ,86)  DISPF 

86  FORMAT! ’FIRST  SELF-SCRLINO  FAILEDI  DISPF  =  ’.1P015.6.  ’  IN”) 
YRRNOE=(Y( ISTOP ) -YD )/( 92 .5/24«D-FtORT ( NCYCLE-1 ) ) 

YMRX=YRRNOE+YtlJN 
YSCRLE=24 .O/YRRNOE 
TOLERrO.O 
C 

C  STACK  DISPLACEMENT  CYCLES  AND  ASSIGN  THE  ROTARY  POT  SIGNAL  TO 
C  THE  VECTOR  RPOT(I) 

C 

9D  TEST=YRRNGE/2.D 
YIM1=Y( ISTRRT) 

RP0T( ISTRRT )=Y( ISTRRT) 

YI=Y( ISTRRT) 

RDDN=O.D 
ISP1=ISTRRT+1 
Y( ISTRRT  )=D  -D 
DO  5  I=ISP1. ISTOP 
RPOT( I  )=Y( I ) 

I F C Y ( I ) .GT .YMAX+TOLER)  GO  TO  101 

IFCYIM1-YCI) -GT.TEST  .RND .  YI.LT.YLIM)  RODN=RDDN+YRRNGE 
YIMlsYCI) 

Y( I  )=Y( I  )+RDDN-YO 
6  YI=Y( I ) 

DO  3  IrlSTRRT. ISTOP 

TM( I )=TIMED+FLORT( I-ISTRRT  JwDELTRT 

1 1  =  1  —  ISTRRT+-1 

DRTRC 1 1 . 1 )=Th( I ) 

DO  3  J=2.NDEG 

DRTRC 1 1 . J )=DRTR( II . J-l )*DRTR( 1 1 . 1 ) 

3  CONTINUE 
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c 

C  ESTIMRTE  NOISE  VARIRNCE  OF  THE  STACKED  SIGNAL  VIA  POLYNOMIAL 
C  REGRESSION 
C 

00  7  I=ISTRRT , ISTOP 
I I=I-ISTART+1 
7  OATHI 1 1 .NDEGP1 )=Y( I ) 

CRLL  MLR( ORTA .X . 1200 .7 .NORTH .NDEGP1 . 1 . .TRUE .  .S .RSQ ) 

VARERR=S 

SD2=SQRT ( VHRERR ) 

SUM=0.0 

SSQ=0.0 

00  10  I  =  ISTRRT » 1 1 
SUM=SUM+Y(  I ) 

10  SSQ=SSQ+Y ( I )**2 

YBAR! LSTART )=SIM/FMRVG 
S2( LSTHRT )=SSQ/FMHVG-YaRR(  LSTART )**2 
C 

C  WRITE  HEADINGS 

C 

C  WRITE! 6 .14) 

14  FORMAT! 1H1 . ’RESULTS  OF  TEST  FOR  OUTLIERS  —  ORIGINAL  VALUES  —  REP 

♦LACEMEWTS ’ ./1H0 .2X» ’ INOEX ’ .T12 . ’POT  VALUE  ( MV )’ ,T29 . ’MOVING  AVG ’ , 
+T46 . ’STO  0EV’.TS8.’REPL  VRL  ( MV )’ ,T76 . ’MOVING  HVG’.T89. 

«’REPL  STD  OEV ’ ) 

C 

C  RECURSIVELY  CALCULATE  MOVING  AVERAGES  ANO  MOVING  VARIANCES  RNO 

C  PURGE  OATH  OF  OUTLIERS.  IF  THE  HIBLER-FILTER  OPTION  IS  SELECTEO. 

C  THE  MOVING  AVERAGE  IS  REPLRCEO  BY  A  LON-PASS  FILTEREO  OUTPUT. 

C 

IF! IHIB.EQ.l )  GO  TO  107 

N0UT=0 

SERR=0 .0 

00  20  L=LSTRRT .LSTOP 

LPRINT=L-ISTART+1 

I0UT=0 

I I=L-ISTHRT+1 
YESTD=X! 1 ) 

00  9  J=2  .NDEGP1 

9  YESTO=YESTO+X( J)*ORTH(II rJ-1) 

YEST(L)=YESTO*YSCRLE 
SERR=SERR+( YESTO-Y ( L ) )**2 
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S01=SQRT(SERR/FL0AT( L-LSTART+1 ) ) 

IF(ABS(Y!L)-YBAR(L) ) .LE.S0N0*SD2)  GO  TO  25 
IOltt=l 
NOLM=NOUT+l 
C 

C  YCL)  IS  AN  OUTLIER.  REPLACE  WITH  Y(L-l)  ANO  CORRECT  MOVING  MEAN 
C 

H0L01=YCL) 

Y!L)=YESTO 
H0L02=YBAR(L ) 

YBAR(L)=YBRR(L)+( Y(L)-H0L01  J/FMAVG 
C  WRITE! 6 . 16 )  LPRINT .HOLD1 .HOLD2.SDI .YCL) .YBAR(L) »S02 
16  FORMAT! 1H  . 15 ,5X . 1P6G15 .6 ) 

26  YBAR!  L+l  )=YBAR( L )+( Y(L+1+MLAG)-Y( L-MLAG)  J/FMAVG 
IFIIOUT)  26.26.20 
26  CONTINUE 

C  NRITEI6.18)  LPRINT. Y(L).YBAR(L). SOI 

18  FORMAT! 1H  . 15 ,5X . 1P3G15 .6 ) 

20  CONTINUE 
WRITE! 6,23)  NOUT 

23  FORMAT! 1H0. ’NUMBER  OF  OUTLIERS  FOUNO  ANO  REPLACED  IS  ’.15) 

24  CONTINUE 

C  WRITE! 6 . 19 ) 

19  FORMAT! 1H1.T6. ’TIME  (S)’.T19.’M  A  VEL  (F/SJ’J 
LG0=LSTART+1 

LSM1=LST0P-1 
00  21  L=LG0,LSM1 

VMA! L )=YSCALE/24 .O/OELTAT*! YBAR! L+l J-YBAR! L-l ) ) 

C  WRITE! 6,22 )  TM( L ) . VMA! L ) 

22  FORMAT! 1H  .1P2G15.6) 

21  CONTINUE 
C 

C  COPY  OUTLIER-PURGEO  OATA  INTO  OATA! I .N0EGP1 ) 

C 

DO  30  I=ISTART , ISTOP 
Y(  I  )=Y! IJxYSCALE 

3D  OATA! I-ISTART+1 .NOEGPl )=Y( I J/12 .0 

CALL  MLR! OATA .X .  1200 ,7 .NOATA .NOEGPl . 1 . .TRUE . ,S ,RSQ ) 

C 

C  WRITE  HEADINGS  FOR  DISPLACEMENT.  VELOCITY  AND  ACCELERATION 
C 


NRITE! 6 ,38 ) 
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38  FORMAT( 1H1 .T6. 'TIME  (S)'.T19.’I  OF  V  (IN )' .T35 . 'VELOCITY  (F/S)'. 
$T50, 'ACCEL  (F/S/S) *.T6’L, 'SIGNAL  (IN) ' .T82. 'RESIOUAL  (FT)') 
0ISPL=0.00  00 
VELT=0.00  00 
00  42  IrlSTART.ISTOP 
I I=I-ISTART+1 
XT=X( 1 )+X( 2 )wOATA( II. 1) 

VEL=X( 2 ) 

ACC=2.0D  OD*X( 3 ) 

DO  40  J=2.NDEG 
XT=XT+X( J+l )*OATA( 1 1 . J  ) 

FJ=FLOAT( J) 

FJJM1=FL0AT( J*( J-l )  ) 

VEL=VEL+FJ»X( J+l )*OATA( 1 1 .J-l ) 

IF(  J.OE .3 )  ACC=ACC+FJJM1*X(  J+l )*OATA(  1 1 . J-2 ) 

40  CONTINUE 

RESID=XT-DATR( 1 1 .N0E0P1 ) 

0ISPL=DISPL+6  .OD  00»0ELTii(  VELT+VEL ) 

VELT=VEL 
XTP=12 .D»XT 

WRITE( 6.43 )  DATA( II .1) .OISPL.VEL.ACC.XTP.RESIO 
43  FORMAT ( 1H  .1P6015.6) 

C 

C  ASSIGN  VELOCITY  TO  YBAR  AND  ACCELERATION  TO  S2  FOR  PLOTTING 
C 

YBAR(  I  )=VE1_ 

S2( I )=ACC 
42  CONTINUE 
C 

C  PLOT  ORIGINAL  POT  OATA 
C 

N  ( 1 )  =  1 
N(  2 )=NDATR 
WRITEt  6.41 ) 

41  FORMAT( 1H1.T20. 'ORIGINAL  ROTARY  POTENTIOMETER  SIGNAL  (MV)') 

CALL  PPLOT( TM( ISTART ) .RPOTt ISTART ) .N .7 .0 .4 -0 .0. '0' ) 

NRITE( 6 .46 ) 

CALL  CPLOT(  TM(  ISTART ) .RPOTt ISTART ) .N.7 .0 .4 .0 .0 .0 .0 . 1 ) 

XPAGE=O.S 
YPAGE=-0 -5 

CALL  SYMBOL (XP AGE. YPAGE. 0.084. 'TIME  (SEC) ' .0.0. 10) 

XPAGE=-0.£ 


52 


YPRGErO.S 

CRLL  SYMBOL (XPRGE. YPRGE, 0.084. ’OISPL  SIO  (MV)  *.90. 0,14) 
XPRGE=0 .5 
YPRGE=RYL+0 .25 

CRLL  SYMBOL( XPRGE .YPRGE .0 .098 .TITLE .0 .0 .NCH ) 

C 

C  PLOT  UNFILTERED  ORTR 
C 

NRITEC  6.44) 

44  FORMAT! 1 H 1 ,T20. 'UNFILTD  DISPLACEMENT  SIGNAL  (MV)’) 

46  FORMRTC 1H  .T20.’TIME  (SEC)*) 

CRLL  PPLOT( TM( ISTRRT).Y( ISTRRT) ,N. 7. 0.4. 0.0, *0’) 
WRITE(6,46) 

CRLL  CPLOTC  TM( ISTRRT ).Y( ISTRRT ) ,N .7 .0 .4 .0 .0 .0 .0 . 1 ) 
XPROE=0 .5 
YPRGE=-0 .5 

CRLL  SYMBOL (XPRGE.YPROE, 0.084, ’TIME  ( SEC )  ’  .0 .0 . 10 ) 

XPRGE=-O.S 

YPRGE=0 .5 

CRLL  SYMBOL (XPRGE.YPROE. 0.084. ’OISPL  SIO  (  MV ) ’ .90 .0 . 14 ) 

XPRGE=0.5 

YPRGE=RYL+0 .26 

CRLL  SYMBOL( XPRGE .YPRGE .0 .098 ,TI TLE .0 .0 .NCH ) 

C 

C  PLOT  ESTIMATE  OF  VELOCITY 
C 

NRITE( 6,74) 

74  FORHRT( 1H1 .T20, ’ESTIMATE  OF  VELOCITY  (F/S)’) 

N( 2 )=NOATR 

CRLL  PPLOT( TM( ISTRRT ) ,YBRR( ISTRRT ) ,N .7 .0 .4.0 .0 . ’V ’ ) 
NRITE( 6.46) 

CRLL  CPLOT( TM( ISTRRT ) ,YBRR( ISTRRT ) .N .7 .0 .4 .0 .0 .0 .0 . 1 ) 
XPRGE=0  *5 
YPRGE=-0 .5 

CRLL  SYMBOL( XPRGE. YPRGE. 0.084. ’TIME  ( SEC) * .0 .0 . 10 ) 

XPRGE=-0 .5 

YPRGE=0.5 

CRLL  SYMBOL(XPRGE. YPRGE. 0. 084. 'VEL  ( F/S ) ’ .90 .0 .9 ) 

XPRGE=O.S 

YPRGE=RYL+0.2S 

CRLL  SYMBOL( XPRGE .YPRGE .0 .098 .TITLE .0 .0 .NCH ) 

C 
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C  PLOT  ESTIMATE  OF  ACCELERATION 
C 

WRITE! 6.76) 

76  FORMAT! 1A-.T20. ’ESTIMATE  OF  ACCELERATION  (F/S/S)’) 
CALL  PPLOT! TM( ISTART) .S2( ISTART) ,W,7 .0 .4 -0 .0 , ’A ’ ) 
WRITE! 6 .46 ) 

CALL  CPLOT! TM( ISTART),S2( ISTART ) .N ,7 .0 .4 .0 .0 .0.0 . 1 ) 

XPAOE=0  *5 

YPAGE=-0.5 

CALL  SYMBOL! XPAGE.YPAGE. 0.084. ’TIME  (SEC)  ’  .0 .0. 10) 
XPAOE=-0  .5 
YPAGE=0 .5 

CALL  SYMBOL! XPAGE.YPAGE, 0.084. ’ACC  (F/S/S )  ’  .90  ,0 . 1 1  ) 

XPAGE=0 .5 

YPAGE=AYL+0.25 

CALL  SYMBOL! XPAGE .YPAGE .0 .098 .TITLE ,0-0 ,NCH ) 

C 

C  PLOT  THE  MOVTNQ  AVERAGE  VELOCITY 
C 

N(2  )=NDATA-2*( MLAG+2 ) 

WRITE! 6 .78 ) 

78  FORMAT!  1H1  ,T20. ’EST.  OF  M-.  A.  VEtOCITY  ( F/S )  ’ ) 

CALL  PPLOT! TM( LOO J.VMA(LGO) ,N .7.0,4 .0.0, ’V’) 

WRITE! 6 . 46 ) 

CALL.  CPLOT!  TM(  LGO  ).  VMAi  LGO  )  ^N, 7. 0.4. 0.0. 0.0.1) 

XPAGE=O.E 

YPAGE=-0 .5 

CALL  SYMBOL! XPAGE ,  YPAGE  .0 .084.  ’TIME  (SEC) ’ .0.0.10) 
XPAGE=-0 .5 
YPAGE=0 .5 

CALL-  SYMBOL! XPAGE.YPAGE. 0.084,  ’VEL  ( F/S ) ’ ,90 .0 ,9 ) 

XPBGE=0 .5 

YPAGE=AYL+0.25 

CALL-  SYMBOL (XPAGE.YPAGE, 0.098, TITLE. O.O.NCH) 

C 

C  SMOOTH  THE  MOViNG  AVERAGE  VELOCITY 
C 

M=6 

C0P=4G.0*DELTAT 
KMRX=2*M+1 
1 1 1=LG0+M 
I  I2=LSM1-M 
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C  NRITE(6.29)  COP 

29  FORMAT!  1H1  ,T10.  'SMOOTHED  MOVING-AVERAGE  VELOCITY  '  ,/lHO.TlO. 

$  ’CUT-OFF  PERIOO  (S)  =  ’ . 1PG15 .6 ./1HO./1HO. T6 ,  ’TIME  (S)’,T20. 
+  'FILT  VEL  ( F/S )  ’  ,T36 . *OISPL  (IN)’) 

00  31  11=111,112 
SUf1=0.0 

00  33  K=1 ,KNflX 

33  SUM=SUf1+C0F2(  K)*VMA(  I  I+K-MLAG-1 ) 

VNBS(  1 1  )=SUt1 

C  NRITE(6,32)  TM.(  1 1 ) . VMAS( II ) ,Y( 1 1 ) 

31  CONTINUE 

32  FORMAJdH  rlP3G15.6) 

C 

C  PLOT  THE  FILTEREO  MOVING-AVERAGE  VELOCITY 
C 

N(2  )=N(2 )-KMAX 
NRITE(6,79) 

79  FORMATdHl  .T20. ’EST.  OF  FILTERED  M.  fl.  VELOCITY  (F/S)’) 

CALL  PPLOT( TM( III) , VhflS( III), N, 7. 0,4. 0,0, ’V’) 

NRITE( 6.46) 

CALL  CPLOT(  TM  ( III) . VMAS(  III), N, 7  *0,4-. 0,0, 0,0,1) 

XPAGE=0 .5 
YPAOE=-0 .5 

CALL  SYMB0L(XPAGE.YPAGE.0.084-.  ’TIME  (SEC)  ’  ,0  .0 . 10 ) 

XPAGE=-0 .5 
YPAGE=O.S 

CALL  S  YMBOL(XPAGE.  YPAGE.  0.084-. ’VEL  ( F/S )  ’  ,90  .0 .9 ) 

XPAGE=0.5 
YPAGE=AYL+0 .25 

CALL  SYMBOL! XPAGE .YPAGE .0 .098 .TITLE ,0 .0 ,NCH ) 

C 

C  PLOT  VEtOCITY  VERSUS  DISPLACEMENT 
C 

NRITE( 6.83) 

83  FORMAT! 1H4  ,T20.  'CROSSPLOT  OF  FILT  VEL  (F/S)  VS  OISPb.(IN)’) 
CALL  PPLOT( YEST( III) ,VMAS( II 1 ) .N.T.0 .4-0 .0 . '\T  ) 

NRITE(6 ,48 ) 

48  FORMAT! 1H  ,T20.’0ISPL  (IN)’) 

CALL-  CPLOT(  YEST(  III) ,  VMAS(  III). N, 7-. 0.4. 0.0, 0.0,1) 

XPAOE=0 .5 
YPAGE=-0 .5 

CALL  SYMBOL( XPAOE .YPAGE .0 .084.  ’DISPL  ( IN)  ’  .0  *0 . 10 ) 
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XPAGEr-0.5 

YPRGE=0.5 

CALL  SYMBOL!  XPRGE  .YPRGE  .0 .084^.  'VEL  (  F/S )  ’  .90  .0  .9 ) 
XPAGE=0 .5 
YPAGE=AYL+0 .25 

CALL  SYMB0L(XPAGE,YPAGE.0.098.TITCE.0.0.NCH) 

C 

C  CLOSE  COMO  FILE  ANO  COPY  TO  OUTPUT 
C 

CRLtr  0$R007(  2  .  * - TERMINAL  DIALOGUE  — M2.0) 

CALL  COMO**! i 000050 .COMOFL .32 .0 .COOE ) 

IF(COOE.OT.O)  GO  TO  103 

CALL  OPENtA!  AtREAO .COHOEL .32 .3 ) 

91  CALL  I$R007( 3 .BUF .40 . $92 ) 

WRITE! 6 . 93 )  BUF 

93  FORMAT!  1H  .40R2) 

GO  TO  91 

92  CALL  CLOSER! 3) 

CALL  OELEtA! COMOFL. 32) 

C 

C  CLOSE  INPUT  ANO  OUTPUT  FILES 

C 

CALL  CLOSER! 1 ) 

CALL  CLOSER! 2) 

C 

C  SPOOL  OUTPUT 

C 

CALL  SPOOL$( 1 . TN. 32. INFO. ORTA. 7168. ICOOE) 

IF! ICOOE.GT.O)  GO  TO  60 
CALL  TNOU!  *  ’.1) 

CALL  TNOUA! ’YOUR  SPOOL  FILE  IS  ’.19) 

CALL  TNOU! INF0(8) .6) 

104  CONTINUE 

C  106  CALL  C*M13( -4.0.NU.0) 

106  CALL  CLOSER! 4) 

CALL  PLOT! 8 .0 ,0 ^0 .999 ) 

CALL  EXIT 

105-  CALL-  TNOU!  ’UNEXPECTED  EOFM  ’ .  16  ) 

GO  TO  106 

107  CONTINUE 
C 

C  FILTER  THE  STACKEO  INPUT  SIGNAL  ANO  PLACE  IN  YELAR(I) 
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c 

CALL  hf ILTR( FC .DELTAT. MLAG. NORTH , ISTRRT .Y.YBAR) 

00  TO  24 
C 

C  PRINT  ERROR  MESSAGE 
C 

60  CRLL  ERRPR$( K$NRTN> ICOOE .0.0. ’SPOOLS ’  ,6 ) 

CALL  EXIT 

101  WRITE!  1.102)  YMAX„YtlIN.YSCALE.I.Y(I) 

102  FORMAT!  ’INCONSISTENT  INPUT  ORTA  * .//  ’YMAX  =  ’ . 1PG1S .6 ./ ’YMIN  =  ’. 
$1PG15*6 ./ ’ YSCREE  =  ’  rIPGIS  .6 ./ *Y(  M3,  * )  =  ’.1PG15.6) 

CALL  TNOUR!  ’ERT  ' 

CRbtT  EXIT 

103  CALL  ERRPRt! KtNRTN .CODE. 0.0 . ’COMO$$ ’  ,6 ) 

END 

SUBROUTINE  C*AD07(NU) 

1  CRLL  I$AD07( NU . IBUF .1.0) 

IF(ANO(  IBUEL,:  177400  )-:000400)  1  .2.1 

2  RETURN 
END 

SUBROUTINE  HFILTR! FC .OELTAT .MLAG.NDATA , ISTART.Y ,Z) 

C 

C  HIBLER’S  NON-RECURSIVE.  GENERAL-PURPOSE.  LOW-PASS  FILTER 

C 

EXTERNAL  $ 

DIMENSION'  C( 101).Y(1).Z(1) 

C 

DATA  PI/3. 1415926636/ 

C 

C  INPUTS: 

C  FC  - 

C  DELTAT  - 

C  MLAG 

C  NORTH  - 

C  ISTRRT  - 

C  Y  - 

C  OUTPUTS: 

C  Z  - 

C 
C 
C 
C 


CUTOFF  FREQUENCY  (HZ) 

TIME  STEP  (SEC) 

MAX  NUMBER,  OF  LAGS  USEO  IN  THE  SYMMETRIC  FILTER 
NUMBER  OF'  DATA  POINTS  IN  THE  (Y)  INPUT 
POSITION  OF  EIRST  OATH  POINT  IN  INPUT  ARRAY 
INPUT  ARRAY 

OUTPUT  ARRAY 

NOTE:  LENGTH  OF  OUTPUT  VECTOR  HAS  BEEN  SHORTENED  BY 
MLAG  POINTS  AT  THE  START  AND  BY  MLAG  POINTS  AT  THE  END 
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EMLAG=FtOAT(  MLAG) 

NC0EFS=2*Mt AG+ 1 
LSTART=ISTART+MLAG 
LSTOP=ISTART +NDATA-1-MLAG 
IF(LSTOP.LE.LSTART)  GO  TO  104 

Nl=2 .0»FMLAG*FC*0ELTAT+0 .001 
IF(Nl.LT.l)  GO  TO  106 
C 

C  CALCULATE  EILTER  COEFS 
C 

DO  10  1=1 .NCOEFS 
N=IA8S( I —ML AO— 1 ) 

EN=ELOAT(  NT 
JMAX=MLAG-1 
SUt1=0.0 

DO  20  Jrl.JMAX 
FJ=EJLOAT(  JT 

SUt1=SUM+H.(  J.N1  )*COS(  PI*FN*FJ/FMLAG) 

20  CONTINUE 

C(  I  )=(  Sl)t1+H(  0  .Ni )/2 .O+HT  MLAG.N1  )/2 .0*C0S(  PI*FNT  J/FtlLAG 
10  CONTINUE 
C 

Cll)rCII)/2.0 

C(  NCOEFS )=C( NCOEFS )/2.0 


C  PERFORM  THE  FILTERING  (CONVOLUTION)  OPERATION 
C 

DO  25  t=LSTART  .LSTOP 
Sl)t1=0.0 

00  26  J=l. NCOEFS 
K=L-MLAG+J-1 
26  Sltt1=SUM+C( J)»Y(  K ) 

26  Z(L)=SUM 
RETURN 

104  WRITEC  1 . 108") 

108“  FORMAT ( 'ERROR  IN  INPUT  TO  SUBROUTINE  HFIETR.  LENGTH  OF  FILTERED  S 
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SERIES  IS  NOT  POSITIVEl ’) 

CALL  EXIT 

106  WRITE! 1,107) 

107  FORMAT!  ’ERROR  IN'  INPUT  TO  SUBROUTINE  HFIhJR .  N1  IS  LESS  THAN  UNIT 
SYI  ’) 

CALL  EXIT 
END 

FUNCTION'  KCI.N1) 

IF(I.OT.Nl)  00  TO  1 

Hal.O 

RETURN 

1  IF! I .EQ.N1  +  1)  H=0 *71 
IF! I .EQ.N1+2)  Hc0.23 
IF! I .CT.N1+2)  H=0 .00 
RETURN 
END 
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Memorandum  for  Record 


SENSITIVITY  OF  INTERIOR  BALLISTICS 
IN  THE  M110A2 
TO 

PROPELLING  CHARGE  TEMPERATURE 


George  Schlenker 
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DRSAR-PEL 


5  December  1979 


MEMORANDUM  FOR  RECORD 

SUBJECT:  Sensitivity  of  Interior  Ballistics  in  the  M110A2  to  Propelling 
Charge  Temperature 


1 .  Reference: 

a.  FONECON  between  Brian  Walters  (DRSAR-HA)  and  George  Schlenker 
(DRSAR-PEL),  29  Nov  79,  subject  as  above. 

b.  Firing  Tables  No.  FT  8-Q-l ,  HQDA,  Jan  76,  title:  Cannon,  8-Inch 
Howitzer,  M201  on  Howitzer,  Heavy,  Self-Propelled,  8-Inch,  M110A1  Firing 
Projectile,  HE,  Ml 06. 

2.  Background 

The  brief  study  reported  in  this  memorandum  was  initiated  in  response 
to  the  Reference  a  phone  conversation.  Although  information  concerning 
the  sensitivity  of  muzzle  velocity  of  the  subject  system  to  propelling 
charge  temperature  is  readily  available,  as  in  Reference  b,  it  did  not 
appear  that  comparable  information  on  sensitivity  of  peak  chamber  pressure 
existed  for  all  charge  zones.  Pressure  information  was  urgently  needed 
to  satisfy  the  requirements  of  government  ammunition  contractors.  This 
MFR  provides  such  estimates  in  a  consistent  manner  for  all  zones  of  the 
Ml  and  M2  propelling  charges  and  for  the  M106  and  M650  projectiles. 

3.  Methodology 

Ballistic  calculations  were  made  using  the  ARRCOM  interior  ballistics 
simulation  for  the  Ml  10/Ml 10A2  howitzer  system.  This  computer  model  had 
been  developed  during  CY  79  to  analyse  the  performance  of  the  howitzer 
with  various  projectiles  fired  from  both  well  seated  positions  and  from 
fallback  positions,  i.e.,  with  the  projectile  resting  upon  the  charge. 

Only  well  seated  projectiles  were  simulated  here.  It  was  anticipated  that 
configurational  differences  in  the  M106  and  M650  projectiles  would  produce 
different  ballistic  sensitivity  to  charge  temperature.  Consequently,  both 
projectiles  were  examined  using  the  Ml  (green  bag)  propelling  charge. 
Additionally,  the  M106  projectile  was  examined  with  the  M2  (white  bag) 
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DRSAR-PEL  5  December  1979 

SUBJECT:  Sensitivity  of  Interior  Ballistics  in  the  M110A2  to  Propelling 
Charge  Temperature 

propelling  charge*.  Peak  chamber  pressure  and  muzzle  velocity  were  obtain¬ 
ed  in  the  M201  cannon  using  these  charges  at  charge  temperatures  of  0,  70, 
and  145  deg  F.  These  results  are  summarized  in  Table  1. 

4.  The  effect  of  propelling  charge  temperature  on  model  parameters  is 
treated  in  the  following  way.  Heat  added  (or  removed)  from  the  charge 
by  raising  (or  lowering)  the  temperature  relative  to  ambient,  70  deg  F, 
is  considered  to  appropriately  raise  (or  lower)  the  flame  temperature 
by  requiring  conservation  of  energy.**  Additionally,  the  effect  of  heat 
addition  (or  removal)  increases  (decreases)  the  linear  burning  rate  of 
the  solid  propellant.  This  phenomenon  can  be  treated  (over  a  restricted 
range  of  temperature)  by  making  the  linear  burning  rate  at  1  ksi,  e»  a 
linear  function  of  temperature.  Thus, 

S(T)  =  e0  ♦  If  (T  -  T0)  . 

where  e(T)  is  the  burning  rate  coefficient  function  of  temperature  T,  and 
and  where  is  the  nominal  value  at  temperature  T  .  The  partial  deriva¬ 
tive,  30/aT,  is  considered  constant  over  the  domain  of  e(T).  However,  it 
is  noted  that  experimental  evidence  indicates  that  3 e/ 3T  decreases  with 
decreasing  temperature.  Comparison  of  the  calculated  with  experimental 
(Reference  b)  decrease  In  muzzle  velocity  with  change  in  charge  tempera- 
ture  from  70  to  0  deg  F  (Table  2)  shows  a  somewhat  larger  calculated 
decrease.  By  inference,  3 3/ 9T  must  decrease  with  decreasing  temperature. 
Unfortunately,  the  form  of  3g/3T  is  not  provided  from  propellant  data 
available  in  the  SPIA  manual.  Therefore,  the  accuracy  of  ballistic  esti¬ 
mates  suffers  for  charge  temperatures  less  than  about  0  deg  F. 


*  Thermochemical  parameters,  propellant  grain  dimensions,  and  charge 
weights  specific  to  certain  charge  lots  were  used  in  the  simulations. 
These  lots  were:  for  the  Ml  charge  IND  69797,  and  for  the  M2  charqe 
BAJ  67951. 

**Energy  conservation  requires  that  the  isochoric  adiabatic  flame  temp¬ 
erature,  V  ,  existing  after  changing  the  charge  temperature  by  AT 

be  related  to  the  ambient  flame  temperature  Ty  by  the  relation: 

rt;/(y  -  d  -  rtv/(y  -  n  +  csolid  at  , 

where  R  is  the  gas  constant  and  y  is  the  ratio  of  specific  heats  for 

the  propellant  gas.  The  specific  heat  of  the  solid  propellant,  in 

consistent  units,  is  C  ...  . 

solid 
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DRSAR-PEL 


5  December  1979 


SUBJECT:  Sensitivity  of  Interior  Ballistics  in  the  M110A2  to  Propelling 


Charge  Temperature 


5.  Resul ts 

The  peak  chamber  pressures  and  muzzle  velocities  calculated  for  all 
computer  runs  are  given  in  Table  1.  Average  values  of  the  pressure 
sensitivity  coefficient  ae/aT  and  of  the  velocity  sensitivity  coefficient 
8V0/3T  are  given  here  for  each  charge  zone.  It  is  generally  a  good  approx¬ 
imation  to  take  these  coefficients  as  constants  for  zones  1G  through  5G 
as  can  be  seen  by  the  linearity  of  p  av  and  V  with  T  in  Figures  1  and  2. 

However,  at  the  higher  zones  using  the  M2  charge  there  is  a  noticeable 
departure  from  linearity.  This  phenomenon  has  been  observed  experimentally 
in  other  systems,  and  is  not  just  an  artifact  of  the  simulation.  Conse¬ 
quently,  the  calculated  average  values  of  ap/aT  and  aV  /3T  for  zones  5W, 

6W,  and  7W  should  not  be  used  for  extrapolation. 

6.  There  are  two  conspicuous  aspects  of  the  numerical  results  which 
deserve  attention:  (1)  Ballistic  sensitivity  increases  in  absolute 
value  with  increasing  zone;  (2)  The  multi -perforated  5W  charge  is  sub¬ 
stantially  more  sensitive  than  the  single-perforated  5G.  One  should  also 
notice  that  the  ballistic  sensitivity  is  less  with  the  M650  projectile 
than  with  the  Ml 06  projectile  at  the  same  zone. 

7.  To  provide  some  measure  of  validity  in  the  calculated  sensitivities, 
a  comparison  is  made  between  calculated  and  nominal  values  of  muzzle 
velocity  sensitivity  obtained  from  the  Firing  Tables  (Reference  b).  One 
possible  reason  for  differences  observed  is  simply  due  to  differences 
between  propelling  charge  lots.  Apart  from  this,  the  assumptions  regard¬ 
ing  the  linear  burning  rate,  mentioned  above,  also  contribute  to  error. 


GEORGE  SCHLENKER 
Operations  Research  Analyst 
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TABLE  1.  CALCULATED  SENSITIVITIES  OF  CHAMBER  PRESSURE  AND  MUZZLE  VELOCITY  TO 
PROPELLING  CHARGE  TEMPERATURE  IN  THE  M201  CANNON  WITH  THE  M106  AND  M650  PROJECTILES 
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TABLE  2.  COMPARISON  OF  CALCULATED  AND  EXPERI¬ 
MENTAL  VALUES  OF  VELOCITY  SENSITIVITY 
TO  PROPELLING  CHARGE  TEMPERATURE 
IN  THE  M201  CANNON  WITH  THE  M106  PROJECTILE 


Charge 

Zone 

Temperature 
(deg  F) 

Calc.  A Vrt 
(f/s)  0 

Tabulated* 

AV 

(f/s) 

V  at  70°F 

0  (f/s) 

1G 

0 

-7.0 

-7.5 

833 

145 

7.8 

9.2 

2G 

0 

-8.0 

-8.5 

912 

145 

8.6 

9.8 

3G 

0 

-8.9 

-9.5 

1012 

145 

9.7 

11.2 

4G 

0 

-10.4 

-10.8 

1156 

145 

11.1 

12.0 

5G 

0 

-12.5 

-12.5 

1388 

145 

13.2 

14.4 

5W 

0 

-27.4 

-23.3 

1464 

145 

27.3 

26. a 

6W 

0 

-32.4 

-27.2 

1708 

145 

32.2 

30.3 

7W 

0 

-37.2 

-30.5 

1995 

145 

37.0 

34.1 

*  Firing  Tables  No.  FT  8-Q-l ,  HQDA,  Jan  76,  title:  Cannon, 
8-Inch  Howitzer,  M201  on  Howitzer,  Heavy,  Self-Propelled, 
8- Inch,  M110A1  Firing  Projectile,  HE,  Ml 06. 
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fH-Mt-lOX  10  TO  1  INCH 
10TH  UNf  HfAVY 
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Pff-MJ-tOX  10  TO  I  MCH 
10THUNI  M«AVY 
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INTERIOR  BALLISTICS 
OF  THE 

Ml 06,  M650  AND  M422  PROJECTILES 
IN  THE 

M110A2  HOWITZER 

WHEN  FIRED  FROM  FALLBACK  POSITIONS 


George  Schlenker 
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DRSAR-PEL 


11  December  1979 


MEMORANDUM  FOR  RECORD 

SUBJECT:  Interior  Ballistics  of  the  M106,  M650,  and  M422  Projectiles  in 
the  M110A2  Howitzer  When  Fired  From  Fallback  Positions 


1 .  Reference : 

a.  Conversations  between  LTC  Stiehl,  DRCPM-NUC-M  (RIA),  and  George 
Schlenker,  DRSAR-PEL,  Nov  79,  subject  as  above. 

b.  Proceedings  of  In-Process  Review,  DRCPM-M110E2,  10  Sep  79,  subject: 
MHO  Cannon  Damage  Study. 

c.  Projectile  Damage  Annex  to  TPR,  DRSAR-PEL,  5  Oct  79,  subject: 

TECOM  Project  No.  2-MU-003-1 06-026. 

d.  Firing  Record  No.  P-82733,  TECOM,  APG,  29  Nov  78,  subject:  Malfunc¬ 
tion  Investigation  Projectile,  8-Inch,  M106  in  M201  and  M2A1  Howitzer  Tubes 
(Projectile  Fallback). 

e.  Report  No.  1375,  MTD  Jefferson  PG,  May  77,  title:  Artillery  Ammu¬ 
nition  Master  and  Reference  Calibration  Chart. 

2.  Background 

The  study  reported  in  this  memorandum  grew  out  of  the  Reference  a  conver¬ 
sations.  The  request  to  study  the  interior  ballistics  of  the  M422  projectile 
when  fired  from  fallback  is  motivated  by  concern  over  the  risk  of  catastrophic 
damage  to  the  projectile  (and  weapon)  when  fired  in  this  manner.  Previous 
emphasis  in  fallback  studies  in  the  Mil  0/Ml 10A2  weapon  system  had  focused 
on  damage  to  the  cannon.  Following  an  in-process  review  (IPR)  of  the  Cannon 
Damage  Study  in  Sep  79  (Reference  b),  increased  attention  was  given  the  risk 
of  damaging  the  projectile  when  firing  from  fallback.  This  increased  empha¬ 
sis  is  illustrated  by  a  proposed  test  program  request  (TPR)  (Reference  c), 
suggested  by  the  M110A2  system  manager,  for  investigating  the  frequency  and 
type  of  damage  incurred  by  each  of  the  following  types  of  projectiles  fired 
from  fallback:  Ml 06,  M650  and  M509.  Although  the  M422  was  omitted  from 
this  TPR,  the  present  study  is  a  first  step  in  assessing  the  risk  of  damage 
to  this  projectile. 
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the  M110A2  Howitzer  When  Fired  From  Fallback  Positions 

3.  Present  evidence  for  the  kind  of  damage  to  the  projectile  expected 
when  firing  from  fallback  comes  from  Reference  d.  In  firing  an  inert 

Ml 06  projectile  in  hard*  fallback  at  zone  2,  smear  camera  coverage  record¬ 
ed  a  separated  fuse  and  a  cracked  projectile  ogive.  From  the  foregoing, 
considerable  interest  has  been  generated  in  describing  the  collision 
between  projectile  and  cannon  when  firing  from  fallback.  This  memorandum 
describes  the  relative  velocity  of  projectile  and  tube  at  first  major 
collision  when  the  area  on  and  near  the  front  bourrelet  first  contacts  the 
rifling.  Also  calculated  and  displayed  is  the  peak  axial  force  arising 
from  this  interaction.  For  comparative  purposes  results  are  obtained  for 
the  three  projectile  types:  Ml 06,  M650,  and  M422. 

4.  Methodology 

Ballistic  calculations  were  made  using  the  ARRCOM  interior  ballistics 
simulation  for  the  MHO  and  M110A2  howitzer  systems.  This  computer  model 
had  been  developed  during  CY  79  to  analyse  the  performance  of  the  howitzer 
with  the  Ml 06  projectile  fired  from  both  well  seated  positions  and  from 
fallback  using  the  Ml  and  M2  propelling  charges.  To  accommodate  the  needs 
of  this  study,  changes** were  made  to  the  program  to  facilitate  the  analysis 
of  other  projectiles  and  other  propelling  charges,  including  the  multi¬ 
grain  M80  charge. 

5.  Phenomena  simulated  in  this  program  include:  (1)  the  burning  of  black 
powder  in  the  base  pad  igniter  with  accompanying  change  in  black  powder 
grain  size  and  unburnt  mass;  (2)  heat  transfer  to  and  ignition  of  the  main 
propelling  charge;  (3)  heat  loss  to  the  projectile  and  cannon  with  conduction 
of  heat  within  the  materials;  (4)  mixing  of  the  gas  components  --  air,  igniter 
gas,  and  charge  gas  --  producing  the  thermochemistry  of  the  mixture;  (5)  gas 
mass  loss  past  the  obturator  (blowby)  throughout  the  interior  ballistic 
cycle;  (6)  burning  of  the  propellant  grain  components  with  associated  grain 
dimensional  change;  (7)  pressurization  of  the  area  behind  the  obturator; 

(8)  development  of  forward  pressure  resisting  projectile  motion  due  to  blow¬ 
by  and  compression  of  the  gas  in  front  of  the  projectile;  (9)  axial  accele¬ 
ration  of  the  projectile;  (10)  acceleration  of  the  unburnt  charge;  (11) 


*  When  fallback  occurs  abruptly  from  a  nearly  seated  position  at  a  high 
quadrant  elevation,  viz,  1150  mils. 

**  Projectile  and  charge  data  were  formerly  represented  in  DATA  statements. 

With  the  requirement  to  simulate  several  other  projectile  and  propelling 
charge  types,  the  program  was  restructured  to  transmit  projectile  and  charge 
data  through  COMMON  statements  from  subroutines  to  the  main  program.  Addi¬ 
tional  program  changes  were  required  to  accommodate  the  use  of  charges 
having  a  mixture  of  grains  from  different  propellant  lots  (having  different 
thermochemistry  ),  as  is  the  case  for  the  M80  propelling  charge. 
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rearward  acceleration  of  the  recoiling  parts;  (12)  development  of  axial 
resisting  forces  on  the  projectile  generated  by  interaction  with  the  gun 
tube.  The  pitch  and  transverse  dynamics  of  the  projectile  are,  specifi¬ 
cally,  not  simulated. 


6.  Due  to  its  extraordinary  nature  in  interior  ballistics,  the  force  of 
interaction  between  projectile  and  tube  during  collision  merits  delineation 
here.  Empirical  evidence  regarding  damage  to  the  lands  of  the  M201  cannon 
indicates  that  the  M106  projectile,  typically,  plastically  deforms  the 
driving  surface  of  the  lands  over  a  printed  area,  A  ,  ,  which  subsequently 

fails  mechanically.  (A  -|n  is  approximately  2  ( i n ^ ) . )  This  collision  area 

experiences  the  material  yield  stress,  oyld>*  during  the  interaction.  Thus, 

an  empirical  estimate  of  the  peak  interfacial  lands  loading  experienced  in 

the  M201  cannon  -  Ml 06  projectile  interaction  is  just 

F  A  ,  cj  *.  i  . 
n  cln  yld 

The  axial  force  associated  with  F  is 

n 

Fzo  =  Fn  (sin  °  +  y  cos 

with  y  the  coefficient  of  friction  and  with  a  the  helix  angle  of  the  rifling 
of  twist  T,  given  by 
a  =  tan”^  (tt/T)  , 

neglecting  slip  with  respect  to  the  rifling.  However,  this  force  can  be 

expected  to  be  a  function  of  the  relative  velocity,  z,  of  the  colliding 

bodies.  To  account  for  this,  the  model  assumes  that  F  occurs  at  an 

zo 

axial  collision  speed  zQ  **  (which  can  be  estimated)  and  that  the  maximum 
axial  force  due  to  an  initial  collision  at  speed  z  is  a  linear  function 
of  z.  Thus, 


Fz  max  Fzo  z^zo  ' 

If  the  projectile  were  being  torqued  to  follow  the  rifling,  an  axial  resist¬ 
ing  force  would  be  generated  F  , 

2  zi 

p  =  2ttp  (sin  a  +  p  COS  a)  z  M 
zl  T  (cos  a  -  y  sin  a)  g~  » 


- 5— - 

ayld  =  Cpst )  >  approximately. 

zQ  =  540  (f/s)  ,  as  given  in  this  program. 
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where 

p  is  the  projectile  radius  of  gyration  in  calibers, 

*z*  the  relative  axial  acceleration,  M  is  the  projectile  weight,  and  g 
is  the  gravitational  constant.  If  the  projectile  angular  velocity  exceeds 
that  required  to  follow  the  rifling  during  collision,  the  axial  force  due 
to  collision  is  assigned  F^ 

7.  Data 


The  author's  previous  experience  with  interior  ballistic  simulation 
supports  the  claim  that  calculated  ballistics  will  be  in  better  agreement 
with  experimental  firings  if  the  specific  chemical  composition  of  a  pro¬ 
pellant  lot  --  as  opposed  to  nominal  composition  --  is  used  to  calculate 
the  thermochemical  parameters  of  the  combustion  products  and  these  para¬ 
meters  are  used  in  the  simulation.  Experience  also  indicates  that  the 
distribution  of  grain  size  actually  measured  from  samples  of  the  propellant 
lot  should  be  used  for  best  results  rather  than  using  a  fictitious  single 
average  grain  geometry.  (The  computer  simulation  actually  uses  nine  grain 
size  classes  obtained  by  treating  grain  dimensions  as  independent  gaussian 
random  variables.)  Lots  of  Ml  and  M2  propelling  charges  have  been  treated 
in  the  manner  described  for  previous  studies. 

8.  For  this  study,  data  for  the  two  propellant  lots  used  in  the  1976  lots 
of  M80  charges  ( PA-76E001 A001 / 2 )  were  used  to  calculate*  the  thermochemi¬ 
cal  parameters  which  characterize  each  propellant  lot.  Results  of  these 
calculations  are  given  in  Table  1. 

9.  Data  on  the  M80  propelling  charge  were, parenthetical ly,  difficult  to 
obtain  and  contained  inconsistencies.  For  example,  the  data  on  pages  48 
through  51  of  Reference  e  indicates  that  the  zone  1  charge  contains  single- 
perforated  propellant  with  an  0.0375  (in)  web,  whereas  information  from 
Radford  Arsenal,  the  propellant  producer,  indicates  that  the  correct  des¬ 
cription  should  have  been  seven-perforated  propellant  with  a  nominal 
0.0360  -  0.0365  (in)  web.  Additionally,  the  web  of  the  MP  zone  2  and 

zone  3  increments  given  in  Reference  e  is  incorrect.  The  value  given 
there  is  0.055  (in)  whereas  the  propellant  lot  acceptance  data  sheets 
show  an  average  web  of  0.068  (in).  Charge  weight  increments  provided  by 
the  PM-NUC  were  also  incorrect.  The  correct  values,  found  in  Table  1, 
were  obtained  from  propellant  acceptance  sheets  and  confirmed  by  lot  firing 


*  Calculations  use  the  Hirschfelder-Sherman  method  with  chemical  constants 
supplied  by  the  Naval  Ordnance  Station,  Indian  Head. 
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records.  Since  the  zones  of  the  M80  charge  consist  of  a  mixture  of  pro¬ 
pellant  lots,  values  of  the  thermochemical  parameters  are  adjusted  for  each 
zone  by  taking  a  mass-weighted  average  of  the  lot-specific  parameters. 

These  results  are  found  in  Table  2. 

10.  During  fallback,  the  projectile  compresses  the  somewhat  loosely  bagged 
propelling  charge.  Due  to  the  great  sensitivity  of  fallback  interior 
ballistics  to  initial  position  of  the  base  of  the  projectile,  it  is  impor¬ 
tant  to  properly  estimate  the  extent  of  propelling  charge  compression.* 

The  degree  of  compression  depends  upon  the  manner  in  which  the  projectile 
falls  back.  If  fallback  is  abrupt  and  occurs  at  a  high  quadrant  elevation 
(QE),  compression  of  the  charge  is  much  greater  than  if  the  projectile 
gradually  slides  back  at  a  low  QE.  The  former  condition  is  referred  to  as 
hard  fallback  and  is  the  condition  simulated  here. 

11.  In  hard  fallback  on  the  Ml  propelling  charge,  the  stresses  induced  in 
the  bag  can  rupture  bag  seams.  However,  when  rupture  does  not  occur,  the 
resulting  charge  deformation  appears  to  occur  at  nearly  constant  charge 
volume,  with  charge  diameter  increasing  to  preserve  volume.  Based  upon 
experiments  with  compression  in  the  Ml  charge,  axial  compression  of  the 
M80  charge  at  constant  charge  volume  is  expected  to  be  limited  by  the 
maximum  diameter  permitted  by  the  chamber.  The  uncompressed  and  compress¬ 
ed  dimensions  of  the  M80  charge  assumed  for  this  study  are  given  in  Table 
3A.  Table  3B  displays  the  empirical  basis  for  the  assumption  of  constant 
charge  volume.  Note  that  the  different  studies  of  charge  compression  have 
yielded  quite  different  results.  This  indicates  that  a  substantial  inherent 
uncertainty  exists  regarding  the  initial  position  of  the  base  of  projectile 
when  fallback  occurs.  In  the  M201  cannon  the  presence  of  standoff  lugs  on 
the  breechface  presents  an  additional  complication.  The  one-inch  deep  lugs 
can  be  expected  to  partially  penetrate  the  charge  during  fallback  compres¬ 
sion.  A  penetration  of  0.5  inch  is  assumed  here. 

12.  The  parameters  which  characterize  the  projectile  in  the  interior  ballis¬ 
tics  program  are  shown  in  Table  4.  Of  particular  significance  for  fallback 
is  the  maximum  diameter  of  the  obturator  and  its  axial  position.  The  leak¬ 
age  of  combustion  gas  is  controlled  by  the  minimum  annular  cross  section 
over  the  projectile  between  the  cannon  internal  profile  and  the  projectile. 


*  The  difference  in  ballistics  associated  with  a  highly  compressed  as 
opposed  to  a  moderately  compressed  charge  is  marked.  In  the  YPG  fall¬ 
back  tests  reported  in  Reference  b,  the  two  conditions  were  produced 
experimentally .  With  the  M106  projectile  and  the  M1/Z2  charge,  hard 
fallback  produced  a  peak  pressure  of  18.6  ksi,  whereas  a  peak  pressure 
of  15.0  ksi  occurred  with  a  nearly  uncompressed  charge  at  200  mils  QE. 
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Generally,  this  throat  occurs  at  the  forward  lip  of  the  obturator  max  dia¬ 
meter.  The  projectile  volume  behind  this  position  is  considered  in  calcu¬ 
lating  the  free  volume  for  the  gas.  All  axial  dimensions  pertaining  to 
the  gun  projectile  interface  --  ZSEAT,  ZCLN  --  are  referenced  to  the  refer¬ 
ence  end  of  the  gun  tube  on  cannon  drawings. 

13.  Results 


An  abbreviated  set  of  ballistic  output  parameters,  calculated  for 
both  fallback  and  well  seated  projectiles  is  displayed  in  Tables  5  and  6. 
For  the  fallback  case  a  relatively  limited  number  of  applicable  experi¬ 
mental  firings  exists.  Peak  chamber  pressure  and  muzzle  velocity  have 
been  the  only  comparable  variables  successfully  measured  in  experimental 
firings.  Attempts  to  measure  collision  accelerations  have  been  unsuccess¬ 
ful.  The  simulated  results  which  can  be  compared  are  in  substantial 
agreement  with  the  pressures  and  velocities  obtained  from  the  Reference  d 
Fallback  Tests  for  zones  2  through  5.  Calculated  velocities  are  generally 
somewhat  in  excess  of  those  measured  (20  to  30  f/s)  but  peak  pressures 
agree  within  the  measurement  error.  As  expected,  there  is  better  agree¬ 
ment  with  experiment  for  the  case  of  well  seated  projectiles. 

14.  The  ballistics  for  well  seated  projectiles  are  offered  here  primarily 
for  comparison  with  those  for  fallback.  For  the  M650  and  M422  projectiles 
there  have  been  no  fallback  firings  to  validate  the  simulated  results. 
However,  these  projectiles  do  not  appear  to  introduce  any  new  phenomena 
with  respect  to  fallback  interior  ballistics.  Consequently,  the  accuracy 
of  results  is  considered  comparable  to  that  with  the  Ml 06  projectile. 

One  should  note  that  the  collision  velocity  at  first  major  collision  may 
not  be  the  best  measure  of  the  severity  of  the  interaction  because  this 
model  does  not  consider  the  effects  of  pitching  and  transverse  motions. 
Using  axial  collision  velocity  as  the  measure  of  severity,  one  can  assert 
that  the  worst  case  for  the  Ml 06  and  for  the  M650  occurs  at  zone  2.  For 
the  M422  projectile  the  worst  case  occurs  at  zone  1  of  the  M80  propelling 
charge.  This  fact  is  due  primarily  to  the  longer  travel  to  first  colli¬ 
sion  at  zone  1  —  13.5  (in)  versus  8.1  (in)  at  zone  2  and  to  the  approxi¬ 
mately  1  (in)  travel  at  start  of  ignition  of  the  M80  propelling  charge. 

15.  Graphical  results  for  the  worst-case  zones  are  included  in  Annex  1. 

A  glossary  of  the  ballistic  variables  displayed  in  Annex  1  is  given  in 
Table  7.  Fallback  simulations  for  the  M106,  M650  and  M422  projectiles  are 
presented  first,  in  that  order  followed  by  the  standard  (well  seated) 
results  for  each  of  these  projectiles  for  the  same  zones.  In  discussing 
the  graphical  results  the  following  points  are  noted.  The  interval  from 
the  time  at  which  the  black  powder  igniter  starts  to  burn  until  the  main 
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charge  is  ignited  is  shorter  for  the  fallback  case  than  for  the  standard 
for  all  projectiles.  The  projectile  has  not  traveled  very  far  at  the 
instant  of  ignition  --  typically,  about  an  inch.  For  all  projectiles 
there  is  a  shorter  Hsetime  to  peak  pressure  for  fallback  than  for  the 
corresponding  standard  case.  The  pressure  and  acceleration  curves  for 
standard  ballistics  are  seen  to  be  quite  smooth  after  shot  start.  How¬ 
ever,  for  the  fallback  case  two  major  disturbances  in  the  smoothness  of 
acceleration  are  evident  following  the  occurrence  of  peak  pressure.  The 
discontinuities  are  due,  first,  to  the  forward  bourrelet  striking  the 
rifling,  which  starts  projectile  spinup,  and,  second,  to  the  rotating 
band  engaging  the  rifling.  There  is  a  significantly  greater  peak  pres¬ 
sure  in  the  fallback  case  than  in  the  standard  case.  For  example,  at 
zone  2  with  the  Ml 06  projectile,  in  hard  fallback  the  peak  pressure  is 
18.8  ksi  versus  10.6  ksi  in  the  standard  case.  With  the  M650  projectile 
using  the  M1/Z2  charge,  fallback  peak  pressure  exceeds  17  ksi  whereas 
standard  peak  pressure  is  about  10  ksi.  The  amplification  of  pressure 
relative  to  standard  due  to  fallback  is  not  as  great  with  the  M422  using 
the  M80/Z1  —  12.0  versus  9.8  ksi.  For  the  Ml 06  and  M422  projectiles, 
the  calculated  muzzle  velocity  is  significantly  higher  (about  50  f/s)  in 
fallback  than  in  the  standard  case.  However,  for  the  M650  projectile  the 
calculated  muzzle  velocity  is  not  greatly  different  in  the  two  cases. 

The  likely  cause  of  the  lower  velocity  for  the  M650  than  for  the  Ml 06 
projectile  at  the  same  zone  in  fallback  is  the  greater  mass  loss  --  1.66 
versus  1.15  lbm  --  due  to  the  smaller  max  obturator  diameter  of  the  M650. 

In  fact,  without  the  large  amount  of  blowby  the  muzzle  velocity  in  fallback 
would  be  much  greater  than  it  is.* 

Summary 

16.  This  study  examines  the  ballistics  of  several  projectiles  -  M106, 

M650,  and  M422  fired  from  both  fallback  and  well  seated  positions  in  the 
M201  cannon.  Motivation  for  this  study  is  a  growing  concern  for  the 
safety  of  these  projectiles  (and  weapon)  when  firing  from  fallback.  The 
ballistic  simulation  produces  estimates  of  the  state  of  the  projectile 
at  first  major  collision  with  the  gun  tube.  The  largest  axial  collision 
velocities  are  produced  with  zone  2  of  the  Ml  propelling  charge  using  Ml 06 
and  M650  projectiles.  These  collision  velocities  are  558  and  486  f/s,  res¬ 
pectively,  the  former  being  sufficient  to  cause  damage  to  the  M106.  For 
the  M422  projectile,  the  worst  case  occurs  with  zone  1  of  the  M80  pro¬ 
pelling  charge  where  the  calculated  collision  velocity  is  389  f/s.  Graphs 


*  Experimental  muzzle  velocities  with  M106  fallback  are  20  to  30  f/s 
greater  than  standard.  Because  of  this  relatively  small  incremental 
velocity  and  a  higher  drag,  the  change  in  range  due  to  fallback  may 
go  unnoticed.  Change  in  range  is  not  a  reliable  indication  of  fallback. 
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of  ballistic  variables  are  presented  in  Annex  1  to  display  the  differences 
between  fallback  and  standard  cases  for  each  of  the  projectile  types. 

s  J-  C  j  v -•*  *■-' 'A--*'  C — ■ 

s'  /• 

1  Incl  GEORGE  SCHLENKER 

Annex  1  Operations  Research  Analyst 


82 


TABLE  1 


COMPARISON  OF  NOMINAL  THERMOCHEMISTRY  OF 
M6  PROPELLANT  WITH  THAT  USED  IN  THE 
M80  PROPELLING  CHARGE 


Propellant  Composition 


Active  Ingredients 

weight  Fraction 

Nominal 

RAD  69535 

RAD  69693 

Nitrocellulose  (13.1 5%N ) 

0.8689 

0.8707 

0.8709 

Dinitrotoluene 

0.0975 

0.0976 

0.0988 

Di butyl  phthalate 

0.0336 

0.0317 

0.0303 

Sum 

1 .0000 

1.0000 

1 .0000 

Additives 

Nominal 

RAD  69535 

RAD  69693 

Diphenyl  amine 

0.0091 

0.0100 

0.0112 

Potassium  sulphate 

0.0099 

0.0109 

0.0110 

Water  (liq) 

0.0050 

0.0080 

0.0044 

Ethyl  alcohol  (resid.) 

0.0090 

0.0021 

0.0134 

Sum 

0.0330 

0.0310 

0.0400 

Thermochemical  Parameter 

Nominal 

RAD  69535 

RAD  69693 

Tv 

(deg  K) 

2549.0 

2597.7 

2517.8 

Y 

1.2544 

1.2524 

1.2552 

n 

(gm  mol/gm) 

0.04405  , 

0.04355  , 

0.04426 

F 

(ft  lbf/lbm) 

0.3122  10° 

0.3146  10° 

0.3099  10 

n 

(in^/lbm) 

30.134 

29.881 

30.232 

M 

(gm/gm  mol) 

22.703 

22.963 

22.593 
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TABLE  2 


VARIATION  OF  THERMOCHEMICAL 
PROPERTIES  WITH  ZONE  IN  THE 
M80  PROPELLING  CHARGE  USING 
PROPELLANT  LOTS  RAD69535  AND  RAD69693 


Parameter  (dimension) 

Zone  1 

Zone  2 

Zone  3 

charge  mass  (Ibm) 
mass  fract.  of  zone  1 

7.456 

15.925 

31.206 

propellant 

1.0000 

0.4682 

0.2389 

Ty  (deg  K) 

2598 

2555 

2537 

Y 

1.2524 

1.2539 

1.2545 

n  (gm  mol/gm) 

0.04355 

0.04393 

0.04409 

F  (106  ft  Ibf  1  Ibm) 

0.3146 

0.3121 

0.3110 

n  (in3/lbm) 

29.881 

30.068 

30.148 

M  (gm/gm  mol ) 

22.963 

22.766 

22.681 

grain  OD  (in) 

0.2065 

0.3745 

grain  PD  (in) 

0.0215 

0.0346 

nom.  web  (in) 

0.0355 

0.0677 

grain  length  (in)  ~ 

0.4800 

0.8610 

std  dev  OD  (%)/(10"Jin) 

3.32/6.8 

1.65/6.2 

std  dev.  len.  (%)/(10~3in) 
a,  burn  rate  expon. 

3,  burn  rate  at  1  ksi 
(in/s) 

9a/9T  (deg  K-1)* 

93/9T  (in/s/deg  K)* 

1.65/7.9 

0.875 

0.235 
-6.81 0-4 
1.010'3 

1 .18/10.2 

*  Inferred  from  strand  burner  data  given  in  SPIA  Manual. 


TABLE  3A 

ESTIMATE  OF  COMPRESSION  IN  THE  M80  PROPELLING  CHARGE 

Dimensions  of  the  M80  Propelling  Charge  Uncompressed  Diameter  7.75  (in) 
Uncompressed  Total  Length  24.25  (in) 


Parameter  Zone  1 

Zone  2 

Zone  3 

Incremental  Length  (in)  5.20 

6.50 

12.55 

Overall  Uncomp.  Length  (in)  5.20 

11.70 

24.25 

Est.  Overall  Comp.  Length  (in)  4.32 

9.73 

20.16 

Initial  Position  in  M201  Cannon  7.6 
(in)  (ZSTART) 

13.0 

23.5 
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TABLE  3B 


ESTIMATES  OF  COMPRESSION  IN  THE 
Ml,  8-INCH  PROPELLING  CHARGE 


Zone 

No. 

Uncompr.  ^ 

Chg.  Length 

Lin) 

Initial  Analytic  Est. 
Compr.  Length 
(in) 

Compressed  Length 
Experimental  Ests 

(2)  (3) 

1 

9.25 

8.25 

7.50 

5.25 

2 

11.00 

10.25 

8.80 

6.25 

3 

13.25 

11.50 

10.75 

11.50 

4 

17.00 

14.50 

13.75 

15.25 

5 

22.00 

19.50 

19.00 

15.00 

(1) 

APG  Firing  Record 

P82733,  Sep  78. 

(2) 

Based  upon  compression  of  lot  IND  69797  with  a 

static  load  of 

150 

lb.  in  a  8.5  (in) 

cylinder.  Reference:  Mario 

Miranda,  YPG, 

Jul  79. 

(3)  Based  upon  the  ARRADCOM  Charge  Damage  Assessment  Tests  using  hard 
fallback  at  1150  mils.  Reference:  Carl  Gardner,  DRDAR-LC,  Aug  79. 


Using  the  uncompressed  length^  and  associated  nominal  max  diameter  of  6.5 
(in)  and  the  experimental  (2)  compressed  length,  the  following  are  the  constant- 
volume  charge  diameters  after  compression: 

Zone  Pi am.  (in) 

1  7.2 

2  7.3 

3  7.2 

4  7.2 

5  7.0 

The  initial  position  of  the  base  of  the  projectile  resting  on  the  propelling 
charge  after  a  hard  fallback  from  a  seated  position  in  the  M201  Cannon  is  esti¬ 
mated  to  be  3.3  (in)  plus  the  compressed  charge  length,  as  measured  from  the 
reference  end  of  the  gun  tube.  This  value  assumes  that  the  standoff  lugs  on 
the  breech  penetrate  to  0.5  (in)  depth  into  the  propelling  charge. 
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TABLE  4 


PROJECTILE  PARAMETERS  USED  IN  THE  INTERIOR 
BALLISTICS  OF  FALLBACK  PROGRAM 


Description  of  Parameter 

Program 

Name 

Ml  06 

Projectile 

M650 

Type 

M422 

Projectile  mass  (Ibm) 

EMP 

200 

199 

242 

Radius  of  gyration  (cal) 

RGRYN 

0.3842 

0.3887 

0.3315 

Dist.  from  base  of  proj.  to 

max  diam.  of  obturator  (in) 

ZBOPRB 

6.49 

3.92 

3.19 

Volume  aft  of  obturator  (in^) 

VBOP 

294.4 

115.2 

141.0 

Max  diam.  of  obturator  (in) 

DIAOBT 

8.28 

8.21 

8.35 

Width  of  band  (in) 

WDBAND 

1.94 

2.04 

3.25 

Dist.  from  obturator  to 

seated  position  (in) 

ZRBSET 

0.00 

2.04 

2.13 

Diameter  of  seat  (in) 

DIASET 

8.28 

8.15 

8.15 

Dist.  from  reference  end  of 
tube  to  base  of  seated  proj. 

(in)  in  M201  cannon 

ZSEAT 

36.30 

38.06 

38.73 

in  M2A2  cannon 

28.26 

30.02 

30.78 

Dist.  from  reference  end  of 
tube  to  base  of  proj.  at 
collision  (in) 

in  M201  cannon 

ZCLN 

28.7 

25.7 

21.1 

in  M2A2  cannon 

20.7 

17.7 

13.2 
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TABLE  5 


SUMMARY  OF  CALCULATED  INTERIOR  BALLISTICS 
IN  THE  M201  CANNON  FOR  SEVERAL  PROJECTILES 
FIRED  FROM  HARD*  FALLBACK 

Simulated  conditions  are  given  below. 


Proj. 

Chg. 

Peak 

Muzzle 

Col  1 ision 

Axial**  Col 

Type 

Type 

Zone 

Press 

Velocity 

Velocity 

Force 

(ksi ) 

(f/s) 

(f/s)  ' 

(k  lbf) 

M106 

Ml 

1 

16.1 

865 

539 

135 

2 

18.8 

959 

558 

140 

3 

21.6 

1066 

552 

138 

4 

26.3 

1218 

523 

131 

5 

31.1 

1428 

275 

no 

M650 

1 

14.5 

815 

473 

119 

2 

17.2 

910 

486 

122 

3 

20.0 

1019 

467 

117 

4 

24.6 

1174 

411 

109 

5 

29.8 

1395 

93 

113 

M422 

M80 

1 

12.0 

936 

389 

109 

2 

19.7 

1392 

329 

110 

3 

36.7 

1995 

0 

19 

*  Abrupt 

drop  from  the  seated 

position  at 

1150  mils 

QE. 

**  The  ratio  of  axial 

to  cross- 

-axial  collision  force 

in  the  M201 

cannon  is 

0.4237 

under  the  assumption 

that  the  projectile  is 

fol lowi ng 

the  rifling. 

Conditions: 

Propelling  charge  temperature  —  70  deg  F. 
Cannon  condition  --  new. 

Ml  charge  lot  no  IND  69797. 

M80/Z1  propellant  lot  no  RAD  69535. 
M80/Z2/Z3  propellant  lot  no  RAD  69693. 
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TABLE  6 


SUMMARY  OF  CALCULATED  INTERIOR  BALLISTICS 
IN  THE  M201  CANNON  FOR  SEVERAL 
PROJECTILES  FIRED  FROM  A  SEATED*  POSITION 

Simulated  conditions  are  given  below. 


Proj. 

Type 

Charge 

Type 

Zone 

Peak 
Press 
(ksi ) 

Muzzle 

Velocity 

(f/s) 

Ml  06 

Ml 

1G 

9.0 

827 

2G 

10.6 

913 

3G 

12.8 

1016 

4G 

16.9 

1164 

5G 

24.5 

1389 

M2 

5W 

14.7 

1477 

6W 

21.2 

1716 

7W 

31.1 

1982 

M650 

Ml 

1G 

8.6 

818 

2G 

10.1 

904 

3G 

12.2 

1005 

4G 

15.8 

1152 

5G 

23.0 

1373 

M422 

M80 

1 

9.8 

874 

2 

14.5 

1279 

3 

29.3 

1880 

*  Projectile  is  assumed  to  have  been  rammed  0.120  (in)  beyond  the 
initial  contact  position. 

Conditions: 

Propelling  charge  temperature  --  70°  F. 

Cannon  condition  --  new. 

Ml  charge  lot  no  IND  69797. 

M2  charge  lot  no  BAJ  67951. 

M80/Z1  propellant  lot  RAD  69535. 

M80/Z2/Z3  propellant  lot  RAD  69693. 
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TABLE  7 


ABBREVIATIONS  USED  IN  PLOTS  OF 
BALLISTIC  SIMULATION  OUTPUT 


Label  Used 

on  Graphs  Detailed  Description 


PRESSURE 

DPDT 

ACCELERATION 
GRAIN  ID  (ENDS) 

GRAIN  OD 

GRAIN  ID  (CNTR) 

CHG  VOLUME 

BURNING  AREA 
MASS  LOSS 
MASS  OF  GAS 
VELOCITY 
INTERNAL  ENERGY 
DISPLACEMENT 

POS  RE  REF 

WALL  TEMPERATURE 
TEMPERATURE 
HEAT  LOSS 
SPIN 


Space-mean  or  chamber  pressure  (psia) 

Pressure  time  derivative  (psi/s ) 

Projectile  axial  acceleration  (g) 

Propellant  gain  inside  diam.  at  ends  of  single 
perf.  grain  (in) 

Propellant  grain  outside  diam.  for  both  SP  and 
MP  grain  (in) 

Propellant  grain  inside  diam.  at  the  center  of 
single  perf.  grain  (in) 

Volume  of  unburnt  propellant  in  propelling 
charge  (in^) 

Area  of  burning  surface  of  propellant  (in^) 

Mass  of  gas  lost  by  escape  past  obturator  (lbm) 
Mass  of  gas  remaining  behind  projectile  (lbm) 
Projectile  velocity  (f/s) 

Internal  energy  of  remaining  gas  (ft-lbf) 
Displacement  of  projectile  relative  to  initial 
position  (in) 

Position  of  the  base  of  the  projectile  relative 
to  the  reference  end  of  the  tube  (in) 
Temperature  at  surface  of  chamber  wall  (deg  K) 
Space-mean  temperature  of  gas  (deg  K) 

Thermal  loss  to  cannon  and  projectile  (ft  lbf) 
Projectile  rotational  frequency  (hz) 
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ANNEX  1 

GRAPHICAL  RESULTS  OF 
INTERIOR  BALLISTICS  FOR 
Ml  06  ,  M  6  5  0  ,  AND  M422  PROJECTILES 
FIRED  FROM  FALLBACK  AND 
WELL  SEATED  POSITIONS  IN 
THE  M  2  0  1  CANNON 
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DRSAR-PEL 


1 3  JAN  1380 


MEMORANDUM  FOR  RECORD 

SUBJECT:  Interior  Ballistics  of  the  M509E1  Projectile  in  the  M110A2 
Howitzer  When  Fired  From  Fallback  Positions 


1 .  Reference: 

a.  MFR,  DRSAR-PEL,  HQ  ARRCOM,  11  Dec  79,  subject:  Interior  Ballistics 
of  the  Ml 06,  M650, and  M422  Projectiles  in  the  M110A2  Howitzer  When  Fired 
From  Fallback  Positions. 

b.  DF,  DRSAR-PE,  HQ  ARRCOM,  20  Nov  79,  subject:  Proposed  Test  Program 
Request  (TPR)  to  Investigate  Projectile  Damage  in  Firing  from  Fallback  in 
the  M110A2  SP  Howitzer,  with  inclosures. 

2.  Background 

This  memorandum  is  a  sequel  to  Reference  a.  Concern  about  the  consequences 
of  firing  projectiles  from  a  fallback  position  in  the  M110A2  howitzer  (M201 
cannon)  extends  to  all  Drojectiles  which  can  be  fired  in  that  system.  Although 
the  original  incentive  for  studying  this  problem  was  to  understand  the  cause 
and  likelihood  of  cannon  damage  --  specifically,  stripped  lands  --  a  more 
recent  emphasis  of  concern  is  to  assess  the  damage  to  the  projectile  fired 
from  a  "hard"  fallback  position.  Hard  fallback  occurs  when  the  nearly  seated 
projectile  abruptly  falls  back  on  the  propelling  charge  during  tube  elevation. 
In  this  condition  the  charge  is  considerably  compressed. 

3.  As  indicated  in  Reference  b,  the  M509  projectile  is  one  of  the  types  of 
projectiles  whose  vulnerability  to  firing  from  fallback  in  the  M201  cannon 
was  to  have  been  examined  experimentally.  The  present  analytic  study  is 
intended  to  support  a  more  extensive  study  of  cannon  and  projectile  damage 
and  to  complement  experimental  work  which  may  proceed  from  the  Reference  b 
proposal . 

4.  Data 


Data  relative  to  the  dimensional  and  inertial  characteristics  of  the 
M509E1  were  received  from  DRDAR-LCU-SS,3  Jan  80,  after  the  Reference  a  study 
had  been  completed.  These  data  reflect  the  most  recent  changes  to  the  aft 
main  body  and  boattail.  Inferences  regarding  the  parameters  used  in  a  computer 
simulation  of  the  M509E1  were  made  from  the  primitive  data.  The  simulation 
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parameters  characterizing  the  M509E1  are  found  in  Table  1. 

5.  Methodology 

The  ballistic  simulation  described  in  Reference  a  was  used  to  calculate 
collision  velocity  and  peak  collision  force  occurring  at  first  major  colli¬ 
sion  between  the  projectile  and  gun  tube  during  a  fallback  trajectory. 

This  event  occurs  when  the  forward  bourrelet  encounters  the  rifling  at  the 
reference  position  ZCLN.  (See  Table  1.)  Another  major  collision  occurs 
when  the  rotating  band  encounters  the  rifling. 

6.  The  dynamic  behavior  of  the  projectile  is  not  described  in  detail  in 
the  ballistics  simulation.  Rather  an  ad  hoc  semiempirical  model  is  used 
to  characterize  the  effect  of  the  collision  on  projectile  (Z-)  axial  and 
spin  motions.  See  Reference  a  for  details.  The  computer  program  reports  a 
variety  of  measurable  variables  such  as  chamber  pressure  and  muzzle  velocity 
as  well  as  many  endogenous  variables  which  would  be  at  best  difficult  to 
measure.  The  simulation  outputs  are  displayed  graphically  for  two  sample 
runs  in  Annex  1.  Ballistics  are  calculated  for  two  conditions:  when  the 
projectile  starts  from  a  hard  fallback  position  and,  second,  when  the  pro¬ 
jectile  is  fully  seated.  The  latter  runs  are  made  for  the  purpose  of  com¬ 
parison  with  the  former. 

7.  Results 


A  brief  summary  of  results  is  provided  in  Table  2  for  the  fallback  case 
and  in  Table  3  for  the  seated  case.  It  is  noted  that  the  largest  collision 
velocity  using  the  M509E1  projectile  occurs  at  zone  1  of  the  Ml  propelling 
charge.  This  contrasts  with  the  situation  using  the  Ml 06  and  M650  projec¬ 
tiles  (Reference  a)  in  which  the  worst-case  collision  velocity  occurs  at 
zone  2.  Further,  differences  are  noted  in  the  magnitude  of  the  collision 
velocity.  For  the  Ml 06  the  largest  velocity  is  558  f/s;  for  the  M650  a 
comparable  value  is  486  f/s;  and  for  the  M509E1  the  largest  collision  velo¬ 
city  is  421  f/s.  Due  to  its  greater  inertia  and  shorter  travel  to  first 
collision,  the  M509E1  suffers  a  less  severe  impact.  Whether  even  this  type 
of  impact  would  damage  the  M509E1  is  not  addressed  here.  It  is  also  noted 
that  the  peak  initial  collision  force  for  the  M509E1  is  nearly  the  same  for 
the  first  four  zones.  This  fact  appears  to  be  due  to  the  combined  effects 
of  projectile  velocity  and  acceleration  during  the  collision,  notwithstanding 
the  fact  that  collision  velocity  declines  with  increasing  zone. 
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8.  One  should  note  from  Tables  2  and  3  that  the  muzzle  velocities  achieved 
with  the  M509E1  are  not  greatly  different  between  the  fallback  and  seated 
cases.  This  is  so  in  spite  of  the  fact  that  the  peak  chamber  pressure  is 
much  larger  with  fallback,  being  almost  twice  as  great  in  fallback  as  stan¬ 
dard  at  zones  1  and  2.  Although  combustion  occurs  at  a  higher  pressure  (and 
with  greater  thermodynamic  efficiency)  in  the  case  of  fallback,  the  losses, 
due  principally  to  blowby,  compensate  for  the  gains  in  this  system. 

9.  Summary 

This  study  examines  the  interior  ballistics  of  the  M509E1  projectile 
fired  in  the  M201  cannon  from  both  fallback  and  well  seated  positions. 
Motivation  for  the  study  is  a  growing  concern  for  the  safety  of  this  pro¬ 
jectile  if  it  were  fired  from  fallback.  The  largest  axial  collision  velo¬ 
city  of  the  M509E1  projectile  is  421  f/s,  which  occurs  at  zone  1  of  the  Ml 
propelling  charge.  Although  this  velocity  is  less  than  the  maximum  of  558 
f/s  produced  with  the  Ml 06  projectile,  a  velocity  of  this  magnitude  may 
still  suffice  to  damage  the  projectile.  Additional  studies  are  required 
to  evaluate  this  possibility.  Graphs  of  ballistic  variables  are  presented 
in  Annex  1  to  display  differences  between  fallback  and  standard  cases  for 
a  typical  firing  zone. 

GEORGE  SCHLENKER 
Operations  Research  Analyst 


1  Incl 
Annex  1 
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TABLE  1 


PROJECTILE  PARAMETERS  USED  IN  THE 
INTERIOR  BALLISTICS  OF  FALLBACK  PROGRAM 


Description  of  Parameter 

Program  Name 

Value  for  M509E1 

Projectile  mass  (Ibm) 

EMP 

207.7 

Polar  radius  of  gyration 

RGRYN 

0.3819 

Distance  from  base  of  proj.  to 

max  diam.  of  obturator  (in) 

ZBOPRB 

4.239 

Proj.  volume  aft  of  obturator  (in 

)  VBOP 

209.7 

Max  diam.  of  obturator  (in) 

DIAOBT 

8.21 

Width  of  band  (in) 

WDBAND 

2.29 

Dist.  from  obturator  to  seated 

position  (in) 

ZRBSET 

1.708 

Diameter  of  seat  (in) 

DIASET 

8.186 

Dist.  from  reference  end  of 
tube  to  base  of  seated  proj. 

(in)  in  M201  cannon 

ZSEAT 

37.74 

in  M2A2  cannon 

29.74 

Dist.  from  reference  end  of 
tube  to  base  of  proj.  at 

collision  (in) 

ZCLN 

in  M201  cannon 

21.52 

in  M2A2  cannon 

13.52 

Plateau  pressure  (psia)  relative 
to  base  of  proj. 

in  M201  cannon 

PPLAT 

450 

in  M2A2  cannon 

400 

Peak  engraving  pressure  (psia) 

PPEAK 

3000 
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TABLE  2 


SUMMARY  OF  CALCULATED  INTERIOR  BALLISTICS 
IN  THE  M201  CANNON  FOR  THE  M509E1  PROJECTILE 
FIRED  FROM  HARD*  FALLBACK 

Simulated  conditions  are  given  below 


Charge 

Type 

Zone 

Peak 

Press 

(ksi) 

Muzzle 

Velocity 

(f/s) 

Collision 

Velocity 

(f/s) 

Axial**  Col 
Force 
(k  lbf) 

Ml 

1 

16.7 

806 

421 

109 

2 

19.6 

900 

412 

109 

3 

22.4 

1009 

357 

109 

4 

27.5 

1164 

232 

111 

5 

32.5 

1380 

0 

24 

*  Abrupt  drop  from  the  seated  position  at  1150  mils  QE. 

The  ratio  of  axial  to  cross-axial  collision  force  in  the  M201  cannon  is 
about  0.4237  under  the  assumption  that  the  projectile  follows  the  riflinq 
Values  of  peak  axial  force  are  given. 

Conditions : 


Propelling  charge  temperature  —  70  deg  F. 

Cannon  condition  *■-  new. 

Ml  charge  lot  no  IND  69797. 

Position  of  base  of  projectile  relative  to  ref.  end  of  tube  after  fallback 
(in)  at  zone: 

1  10.8  (in) 

2  12.2 

3  14.3 

4  17.0 

5  22.3 
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TABLE  3 


SUMMARY  OF  CALCULATED  INTERIOR  BALLISTICS 
IN  THE  M201  CANNON  FOR  THE  M509E1 
PROJECTILE  FIRED  FROM  A  SEATED*  POSITION 

Simulated  conditions  are  given  below. 


Charge 

Type 

Zone 

Peak 

Press  (ksl) 

Muzzle 

Vel  (f/s) 

Ml 

1G 

8.7 

808 

2G 

10.3 

892 

3G 

12.5 

993 

4G 

16.4 

1138 

5G 

23.9 

1358 

M2 

5W 

14.5 

1442 

6W 

20.9 

1676 

7W 

30.9 

1938 

Projectile  Is  assumed  to  have  been  rammed  0.120  (in)  beyond  the  initial 
contact  position, 

Conditions : 


Propelling  charge  temperature 
Cannon  condition  -  new. 

Ml  charge  lot  no  IND  69797. 

M2  charge  lot  no  BAJ  67951. 


70  deg  F. 
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ANNEX  1 

GRAPHICAL  RESULTS  OF 
INTERIOR  BALLISTICS  FOR  THE 
M509E1  PROJECTILE  FIRED  FROM 
FALLBACK  AND  WELL  SEATED 
POSITIONS  IN  THE  M201  CANNON 
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TABLE  A1 


ABBREVIATIONS  USED  IN  PLOTS  OF 
BALLISTIC  SIMULATION  OUTPUT 


Label  on 
Graph 


Description 


PRESSURE 

DPDT 

ACCELERATION 
GRAIN  ID  (ENDS) 

GRAIN  ID  (CNTR) 

GRAIN  OD 
CHG  VOLUME 

BURNING  AREA 
MASS  LOSS 
MASS  OF  GAS 
INTERNAL  ENERGY 
VELOCITY 
DISPLACEMENT 

POS  RE  REF 

WALL  TEMPERATURE 

TEMPERATURE 
TOT  FWD  RESISTANCE 

ENGR  RESISTANCE 

HEAT  LOSS 
FILM  COEFFICIENT 

SPIN 


Space-mean  or  "chamber"  pressure  ( psi a ) 

Pressure  time  derivative  (psi/s) 

Projectile  axial  acceleration  (g) 

Propellant  grain  inside  diam.  at  the  ends  of 
single  perf.  grain  (in) 

Propellant  grain  inside  diam.  at  the  center  of 
single  perf.  grain  (in) 

Propellant  grain  outside  diam.  (in) 

Volume  of  unburnt  propellant  in  propelling  charge 
(In3) 

Area  of  burning  surface  of  propellant  (in^) 

Mass  of  gas  lost  by  escape  past  obturator  (Ibm) 
Mass  of  gas  remaining  behind  projectile  (Ibm) 
Internal  energy  of  remaining  gas  (ft  Ibf) 
Projectile  axial  velocity  (f/s) 

Axial  displacement  of  projectile  relative  to 
initial  position  (in) 

Position  of  the  base  of  the  projectile  relative 
to  the  reference  end  of  the  tube  (in) 

Temperature  at  surface  of  the  chamber  wall 
(deg  K) 

Space-mean  temperature  of  the  gas  (deg  K) 

Total  resistance  to  projectile  forward  motion 
(Ibf) 

Resistance  to  projectile  forward  motion  due  to 
band  engraving  and  friction  (Ibf) 

Thermal  loss  to  cannon  and  projectile  (ft  Ibf) 
Film  coefficient  of  heat  transfer  (ft  Ibf/ 
in2/s/deg  F) 

Projectile  rotational  frequency  (hz) 
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DRSAR-PEL 


2  November  1979 


MEMORANDUM  FOR  RECORD 

SUBJECT:  Approximation  for  the  Standardized  Cumulative  Production  Cost 
Using  Learning  Theory 


1.  Reference: 

a.  Handbook  ALM-63-31 26-H2,  US  Army  Logistic  Management  Center  (ALMC), 
issued  first  by  Directorate  of  Procurement  and  Production,  USAMICOM,  title: 
Alpha  and  Omega  and  the  Experience  Curve. 

b.  CEE  Course  Handout  ALM-63-31 26-H1C,  US  ALMC,  title:  Learning  Curve 
Summary. 

2.  Background 

I  recently  had  the  opportunity  to  attend  a  course  offered  by  ALMC  which 
was  concerned  with  cost  estimation.  During  the  course,  materials  were  pre¬ 
sented  concerning  the  theory  of  learning  or  experience  related  to  product 
manufacturing  cost.  This  MFR  was  stimulated  by  that  course. 

3.  Where  applicable,  learning  theory  is  useful  in  predicting  the  relation 
of  manufacturing  unit  cost  (or  man  hours)  of  a  product  to  the  number  of 
units  produced.  Manufacturing  experience  with  missiles  and  with  tanks 
and  automotive  vehicles  has  shown  that  the  specific  reduction  of  unit  cost 
with  units  produced  is  commodity  dependent.  However,  for  a  variety  of  pro¬ 
ducts  and  over  a  considerable  range  of  production  quantity,  the  expected 
direct  cost  of  manufacturing  the  2n  th  unit  is  proportional  to  the  unit 
cost  for  the  n  th  unit  for  n  an  integer  up  to  some  limit,  where  the  constant 
of  proportionality,  s,  is  called  the  "slope"  of  the  learning  curve  and  is 
frequently  expressed  as  a  percentage.*  Thus,  the  unit  cost  for  the  n  th 
unit,  u(n),  is  given  by 

u(2n)  =  s  u(n) ,  0<S£l  .  (1) 

This  result  implies 

u(n)  =  AnB,  (2) 


Production  experience  with  missiles  indicates  a  79  to  84%  slope  and  with 
automotive  and  tank  vehicles  a  92  to  96%  slope. 
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Cumulative  Production  Cost 


with  A  the  first  unit  cost  and 


B  =  log(s )/log(2)  .  (3) 

From  (2)  the  expected  total  cost  of  manufacturing  a  lot  of  products 
starting  with  the  n1  th  and  proceeding  through  the  n2  th  is  given  by 


total  lot  cost  =  zn^u(n) 
or 

n 


"  A[V2/  -  • 


(4) 

(5) 


The  sum  terms  on  the  right  in  equation  (5)  have  been  named  cumulative 
total  factors  for  the  n2  th  and  n]  -1  th  quantities.  Because  the  term 

Q<n> =  4T/  (6) 

represents  the  cumulative  production  cost  from  the  first  through  the 
N  th  unit  standardized  by  division  by  the  first  unit  cost,  it  is  here  called 
the  standardized  cumulative  cost.  With  this  notation  equation  (5)  becomes 

total  lot  cost/A  =  Q(n2)  -  Q(n,-1)  .  (7) 

4.  Appl ications 


A  variety  of  applications  of  (7)  are  possible.  For  example,  one  may  have 

prior  estimates  of  A  and  B  and  wish  to  calculate  the  total  lot  cost  or  averaqe 
lot  cost  —  J 


total  lot  cost/(n2  -  n]  +  1 ) 

Alternatively,  one  may  have  a  total  lot  cost  and  an  estimate  of  B  from  prior 
experience  and  wish  to  estimate  A.  Finally,  one  may  wish  to  estimate  the 
learning  slope  B  from  values  of  total  lot  cost  and  first  unit  cost.  Addi¬ 
tionally,  the  "algebraic  lot  midpoint",  i.e.,  the  value  of  n  (=K)  such  that 

AKB 


is  the  average  lot  cost,  is  obtained  from  (7)  via  the  equation 


„B  ^n2}  "  ^V1) 

K  -  — - - - 

n2  -  n,  +  1 


AQ 

n2  -  n,  +  1 


(  o\ 

V'j 
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5.  Because  of  the  significance  of  the  standardized  cumulative  cost  and  the 
effort  to  calculate  it  directly  from  (6),  this  quantity  has  been  tabulated 
for  a  large  set  of  the  parameters  B  (or  s)  and  n.  Some  examples  of  tables 
are  given  in  Reference  a.  Without  recourse  to  tables  or  a  computer,  one 
can  calculate  Q(n)  using  a  numerical  approximation. 


6.  Approximations  for  Q(n) 


One  approximations  for  Q(n)  is  obtained  directly 
for  the  algebraic  lot  midpoint,  ,  proposed  by  ALMC 


from  the  approximation 
in  Reference  b: 


_  F  +  L  +  2/FL 

1  4 


(9) 


with 


F  =  first  unit  in  the  lot 
L  =  last  unit  in  the  lot. 

With  F  =  1  and  L  =  n  and  using  the  definition  of  algebraic  lot  mid¬ 
point,  the  ALMC  approximation  for  Q(n)  is 

Q-j  (n)  =  nK.|B 

=  n[(l  +  n  +  2/in)/4]B  .  (10) 

7.  Another  approximation,  proposed  by  the  author,  is  derived  by  treating 
n  as  a  continuous  variable  and  by  solving  the  integral  analog  of  equation 
(6).  This  approximation,  designated  Q2(n),  has  the  form 

Q2(n)  =  (B  +  l)‘1[(n  +  0.5)B+1  -  0.5B+1]  .  (11) 

The  presence  of  the  term  0.5  as  a  correction  to  n  is  motivated  by  the 
need  for  a  symmetric  integral  approximation  to  a  discrete  variable  in  the 
same  spirit  as  the  gaussian  distribution  function  approximates  a  binomial 
distribution. 

8.  Using  the  approximation  Q2(n),  the  standardized  cumulative  cost  of 

manufacturing  a  product  from  the  n.  th  through  the  n?  th  units  is  given 
approximately  by  1  y 

AQ  =  (B  +  l)_1[(n2  +  0.5)B+1  -  (n]  -  0.5) 8+1 \  .  (12) 
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In  the  limit  as  B  approaches  zero,  i.e.,  as  no  learning  occurs, 
aQ  =>  n^  -  n.j  +  1 , 

which  is  exactly  what  would  be  expected.  Further,  the  standard¬ 
ized  average  unit  cost  of  manufacturing  n  units  from  the  n^  th  through 
the  n^  th: 


n. 


n^  +  n  -  1 


15  9lVe"  b*  „  <„,  ♦  0.5)B+1  -  («.  -  0.5)B+1 

lot  std  avg  unit  cost  -  (B  +  T y (n?  .  \  , , - 


(13) 


9.  Numerical  Comparisons  of  Accuracy 

With  the  availability  of  scientific  calculators  the  two  approximations 
given  in  (10)  and  (11)  can  be  easily  carried  out.  Some  numerical  examples 
are  provided  in  Table  1.  The  range  of  slope  values  shown  there  is  repre¬ 
sentative  of  a  variety  of  products.  Although  lot  sizes,  n,  may  exceed  1000 
units,  Table  1  does  not  display  Q(n)  for  larger  values  since  the  theoretical 
or  ideal  learning  may  have  leveled  off  by  that  point.  Note  that  for  values 

A  A 

of  n  less  than  about  20,  and  are  equally  good  approximations.  However, 
for  n  >100,  Q„  is  distinctly  superior  to  the  ALMC  approximation,  Q. .  Note 
also  that  the  error  in  (100)  is  about  2.4%  whereas  that  of  Q2 ( 1 00 )  is  only 
0.02%  for  a  90%  slope.  The  approximation  Q2(n)  improves  as  n  increases, 
while  Q-j  (n )  becomes  a  progressively  poorer  approximation  with  increasing 
n.  The  foregoing  results  are  generally  valid  for  learning  slopes  above  80%. 

10.  Summary  and  Recommendations 

It  has  been  shown  that  application  of  learning  theory  requires  the  use 
of  the  standardized  cumulative  cost  (sometimes  called  the  cumulative  total). 

To  retain  tables  of  this  quantity  is  unnecessary  if  one  has  use  of  a  present¬ 
ly  commonplace,  scientific  calculator.  The  use  of  an  approximation  for  Q(ri) 
seems  justified  since  the  error  of  the  approximation  for  the  best  approxi¬ 
mation  Q2(n)  is  quite  small  relative  to  other  errors.  Since  Q9  is  just  as 
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/V  A 

easy  to  calculate  as  Q, ,  and  since  it  is  generally  more  accurate, 

■  a  L. 

is  recommended.  Further,  Q.|  is  to  be  avoided  for  n  in  excess  of  200 
because  of  a  sizeable  error  of  approximation. 


c  l  j  '  —  ■A-'/i.  <■  -•  «•  L 

y 

1  Incl  GEORGE  SCrfLENKER 

Table  1  Operations  Research  Analyst 
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TABLE  1.  COMPARISON  OF  APPROXIMATIONS  FOR  THE  STANDARDIZED  CUMULATIVE  MFG 
COST  WITH  THE  EXACT  VALUE  FOR  SEVERAL  LEARNING  CURVES 


Percent 

Learning 

Slope 

B* 

Number 

of 

Units 

Cumulative  Cost 
Approx  (1)  A  (2) 

«l  Q2 

Exact  (3) 
Val  ue 

80 

-0.321928 

5 

3.67 

3.76 

3.74 

10 

6.24 

6.34 

6.32 

20 

10.46 

10.51 

10.48 

50 

20.36 

20.15 

20.12 

100 

33.37 

32.68 

32.65 

200 

54.32 

52.75 

52.72 

1000 

165.70 

158.70 

158.67 

85 

-0.234465 

5 

3.99 

4.05 

4.03 

10 

7.09 

7.13 

7.12 

20 

12.48 

12.42 

12.40 

50 

25.99 

25.53 

25.51 

100 

44.96 

43.77 

43.75 

200 

77.40 

74.81 

74.79 

1000 

270.04 

257.94 

257.91 

90 

-0.152003 

5 

4.32 

4.35 

4.34 

10 

8.00 

8.00 

7.99 

20 

14.73 

14.62 

14.61 

50 

32.72 

32.15 

32.14 

100 

59.56 

58.15 

58.14 

200 

108.08 

104.98 

104.96 

1000 

427.95 

412.18 

412.17 

95 

-0.0740009 

5 

4.66 

4.67 

4.66 

10 

8.97 

8.96 

8.95 

20 

17.23 

17.14 

17.13 

50 

40.67 

40.23 

40.22 

100 

77.70 

76.59 

76.58 

200 

148.22 

145.70 

145.69 

1000 

661.53 

647.45 

647.42 

Notes : 

B  =  log10(slope)/log10(2) 

g 

(1)  Cumulative  using  the  ALMC  algebraic  lot  midpoint  is  n[—  ■*  n.+ 

(2)  Q2  =  (B  +  1 )_1[(n  +  0.5)B+1  -  0. 5B+1 ] 

(3)  Q=Ej=1iB 


Incl  1 
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